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Abstract. This course is designed to give a mathematically coherent 
introduction to the classical thory of black holes and also of strings and mem- 
branes (which are like the horizon of a black hole in being examples of physical 
systems based on a dynamically evolving world sheet) giving particular atten- 
tion given to the study of the geometry of their equilibrium states. 



Preface. 

The purpose of this course is to provide a mathematically coherent 
introduction to the classical theory of black holes and also to the related and 
more recently developed topic of the classical theory of relativistic strings and 
membranes for which many of the same techniques are required. The stategy 
of the course will be to concentrate on general results rather than special 
examples, and to distinguish as clearly as possible what has been completely 
proved from what has only been partly established or merely conjectured 
so as to give some idea of the main open problems for future research. The 
discussion is developed on the basis of a chain of key results for which it 
has been possible to provide reasonably complete and self contained math- 
ematical proofs without resort to disproportionate technical complication. 
The level of previous knowledged required corresponds to what is obtainable 
from the relevant sections of a textbook such as that of Misner, Thorne and 
Wheelert^l (whose notation will be used as far as possible) or, in a less ency- 
clopaedaic but more conveniently accessible (and up to date) form, that of 
Waldl^l. There are already several textbooks specifically devoted to various 
aspects of black hole theory [3l[4][5][6] . attention is particularly to be drawn 
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that of Hawking and Ellist^^ for advanced mathematical background read- 
ing, and to that of Novikov and Frolov^^^ for an exceptionally comprehensive 
survey of the published litterature including more than 600 references. 
The organisation of the course is as follows. 

Section 1 provides a brief astrophysical introduction consisting essen- 
tially of a simple explanation'^! of the orders of magnitude that are relevant 
to the conventional idea of the formation of "ordinary" black holes by stellar 
collapse (no such simple and clear picture being available for the more exotic 
phenomenon of the giant black holes that are commonly believed to be lie at 
the heart of active galactive nuclei). 

After this physical introduction, the main part of the course is more 
essentially mathematical in nature following more or less the same lines as 
my previous reviews'^^^t^^^t^^! though with the omission, except for the neces- 
sary references, of certain parts in order to make way for the inclusion of new 
results. Section 2 presents some of the main results of the theory of exactly 
spherical gravitational collapse, which is the only case for which a precise 
dynamical analysis is available. Section 3 gives a brief account of what little 
is known about dynamical formation of black holes in more realistic situa- 
tions where spherical symmetry is broken by effects such as rotation. Section 
4 deals with the theory of stationary rotating black hole equilibrium states in 
the general case for which externally orbiting matter rings may be present. 
Section 5 deals more specifically with the uniqueness theorem that is avail- 
able when no external sources are present. Section 6 concludes the course on 
black hole theory by describing some of the rather miraculous special proper- 
ties of the ensuing Kerr Newman metrics, whose stability is one of the most 
important topics that (for lack of time and space) has not been included in 
this course : for the most complete result, going a long way towards confirm- 
ing that these equilibrium solutions can indeed be considered to be stable, 
the interested reader is referred to the recent work of Whiting'^^l 

Section 7 moves on to present a covariant formulation'^^' of the ba- 
sic mechanical principles of classical brane theory meaning the subject that 
includes the theories of point particles, strings, membranes and continua as 
special cases. Section 8 deals more specificly with the theory of spacially 
isotropic branes, a category that includes all classical string models and in 
particular those representing "superconducting cosmic strings" . 

Finally in a purely mathematical appendix, some of the most impor- 
tant tensorial quantities (which are useful for black hole theory and indispen- 
sible for brane theory) characterising the different kinds of curvature of an 
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imbedding are presented in a readily utilisable form'^^^ that is not yet readily 
available elsewhere. 

1. The astrophysical context of Black Hole formation. 

The study of black holes in general, and of black hole equilib- 
rium states in particular, arises as a natural offshoot of the study of stellar 
equilibrium states whose theoretical foundations were established by workers 
such as Eddington and Chandrasekhar in the years following the elucidation 
of the basic principles of quantum mechanics. In terms of the fundamental 
Plank type unit system that will be used throughout this course (in which 
the speed of light c, Newton's gravitational constant G, the Dirac - Plank 
constant h and the Boltzman constant k are all simuUtaneously set equal 
to unity) the dominant physical mechanisms governing the situation can be 
described^^ll^^l in crude order of magnitude terms (give or take a power of 
ten here or there) in terms of just three particularly important dimensionless 
parameters, namely the masses me and of the electron and the proton, 
and the magitude e of their electric charge, which are expressible as the mod- 
erately small "fine structure" coupling contant ~ 1/137 the considerably 
smaller mass ratio mg/mp ~ 1/1800 and the extremely small gravitational 
coupling constant 10~^^. 

In the low temperature limit, the equilibrium states of small, 
medium, and even moderately large bodies, on scales ranging from single 
molecules through sand grains up to entire planets, are characterisable in 
crude order of magnitude by a typical density p given by 

p e^rrie^mp (1.1) 

which works out (by no means accidentally) to be very roughly of the order 
of unity in "ordinary" units, gm/cm? (which have of course been deliberately 
normalised to give such a result). Taking account of the fact that in all such 
states the mean mass per baryon is given to a very good (within one per cent) 
accuracy by m^, so that p ^ rrtpn where n is the baryon number density, the 
relation (1.1) expresses the condition that the mean separation A ^ n~^/^ 
between baryons will be of the same order as that between the (within a 
factor of two equally numerous) electrons, and therefore of the order of the 
Bohr radius, A ~ l/e^rrie, which is the result that is obtainable from the con- 
sideration that the equilibrium is determined by the balance between Fermi 
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(exclusion principle) repulsion between electrons and electrostatic attraction 
between negatively charged electrons and positively charged ions. 

Although applicable to bodies on scales ranging from that of a 
hydrogen atom to that of the earth, the formula (1.1) loses its validity for 
bodies so large that the long range cumulative effect of the (individually 
very weak) gravitational attraction forces becomes stronger than the effect 
of the electrostatic attraction forces (which of course only act locally because 
of the long range cancellation resulting from overall electric neutrality). For 
a body of mass M, mean density p and hence characteristic mean radius 
R ~ {M/pY^^ resistance to collapse under the influence of gravitational self 
attraction requires a mean central pressure P given according to the well 
known "virial theorem" by 

P^M^V^^^ (1.2) 

which expresses a balance between the typical radial pressure gradient, of 
order P/R, and the gravitational force density, of order pM/R"^. 

The pressure contribution resulting from the application of the 
Fermi exclusion principle to the electrons is of the order of the corresponding 
kinetic energy density, and therefore will be given - in the non relativistic 
limit - roughly by 

p^^ur (1,3) 

me \mpj 

in view of the fact that the mean momentum per electron will just be the 
inverse of the corresponding De Broglie wavelength, which will itself be 
of the same order of magnitude as the mean separation, A ~ n~^^^, where 
n c:^ p/ rUp. So long as the mass M is small compared with a critical value given 
roughly hy M !^ &" jm^ ^ the virial pressure requirement (1.2) is small com- 
pared with the Fermi energy density (1.3) at the "ordinary" matter density 
(1.1) which means that the gravitational compression effect will be unimpor- 
tant, but but beyond this critical mass (which is of the order of that of the 
giant planet Jupiter) the long range gravitational attraction will dominate 
over the short range electrostatic binding so that the corresponding equilib- 
rium states will be of white dwarf type, with the central pressure determined 
by direct equation of (1.2) and (1.3) which means that the characteristic 
mean central density p will be given as a function of the mass M by an order 
of magnitude relation of the form 

p ^ m^vTipM'^ . (1.4) 
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The range of vahdity of the relation (1.4) is of course limited to 
that of the non relativistic degenerate electron gas pressure formula (1.3) 
from which it is derived. When the relevant DeBroglie wavelength A ~ n~^^^ 
becomes short compared with the Compton wavelength A ~ mg, the kinetic 
energy per electron is no longer given by l/A^mg but just by 1/A, so that 
the non relativistic formula (1.3) must then be replaced by the corresponding 
relativistic degenerate gas pressure formula 



(1.5) 

The (by now generally accepted) recognition that the theory of black holes 
must be taken seriously as something whose implications are directly relevant 
and testable in observational astrophysics derives from the startling (1930) 
discovery by Chandrasekhar'^''! that substitution of (1.5) instead of (1.4) in 
the virial equilibrium condition (1.2) does not just give a modified version 
of the functional relation (1.4) for the equilibrium density p as a function of 
the mass M, but instead gives an absolute cut off at a critical mass 

(1.6) 

above which no ordinary cold equilibrium state is possible at all ! 

The existence of this upper mass limit does not of course mean 
that there are no cold equilibrium states beyond the critical density p ~ 
me^mp at which the white dwarf range (1.4) reaches the Chandrasekar limit 
(1.6), since it is also possible to have high density states in which the electrons 
are combined with protons to form neutrons for which the relevant analogue 
of the non relativistic degenerate gas pressure formula (1.3) is 

p^^i^r . (1.7) 

nip \mpj 

However the resulting range of neutron star equilibrium states, with density 
given, by substitution of (1.7) in (1.2), as 

p ^ rUp^M^, (1.8) 

will be cut off by an upper mass limit that is still given'^^^l'^^l by the same crude 
order of magitude formula (1.6) as before, because the relativistic degenerate 
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gas pressure has the same form (1.5) for neutrons as for electrons. More exact 
calculations (whose results are still subject to a considerable uncertainty 
due to the imprecision of our present understanding of the detailed physical 
properties of neutron star matter) indicate that the upper mass limit for 
neutron stars is somewhat larger (though only by a modest factor not much 
in excess of two), than that for white dwarfs : this conclusion is of great 
astrophysical importance, and would appear to have been observationally 
confirmed by the discovery of pulsars (since if the exact upper mass limit 
for neutron stars had turned out to be smaller than that for white dwarfs 
then the formation of neutron stars by gravitational collapse would have been 
rendered virtually impossible). 

As the astronomical community belatedly recognised (after more 
than thirty years of general indifference or incredulity) Chandrasekhar's discovery'^^^ 
made it absolutely necessary to take the possibility of runaway gravitational 
collapse - and ensuing formation of massive or ultramassive black holes - very 
seriously as a phenomenon of potentially crucial relevance to many directly 
observable phenomena. This contrasts with the situation that still applies to 
speculation on the subject of microscopic black holes (for which quantum 
phenomena such as Hawking radiation'^^l are significant) whose relevance to 
anything actually observable remains subject to reasonable doubt. Neverthe- 
less the existence of several categories of observational "black hole candi- 
dates" (of which the most famous prototype example is the galactic X-ray 
source Cygnus X-1) does not yet amount to a firm confirmation of that black 
holes with the properties described in the following sections of this course re- 
ally do exist. The most numerically numerous (and perhaps ultimately most 
atrophysically important) category of observationally detected "candidates" 
is that of nuclei of "active" galaxies, but such (ultramassive) objects are all 
too fuzzy and far away to have been of any use so far from the point of 
view of verification of the basic physical theory. As far as the more conve- 
niently tractable candidates within our galaxy are concerned, the awkward 
fact to be faced is that sixty years after our colleague Chandra's precocious 
and revolutionary theoretical discovery, and more than twenty years after 
the collapse of psychological resistance to the notion of a "black hole" follow- 
ing the coining of the term itself (by John Wheeler) and the (approximately 
simultaneous and no less psychologically significant) experimental discovery 
(by Jocelyne Bell and Tony Hewish) of the pulsars whose identification as 
neutrons stars has long been unquestionable, there are still disappointingly 
few observationally discovered objects that can plausibly be interpreted as 
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"ordinary" (moderate sized) black holes. 

The relative scarcity of black holes in the mass range immediately 
above the Chandrasekhar limit (which is about one and a half times the 
mass of the sun) might at first seem surprising in view of the fact that this 
particular mass range is precisely that of the most numerous subclass of the 
ordinary stars that are visible at night to the naked eye. However from a 
theoretical point of view this apparent paradox can easily be understood as 
follows. 

The non-existence of any material "ground state" configurations, 
i.e. minimum energy (cold, static, absolutely stable) equilibrium states, above 
a critical mass value given in order of magnitude by (1.6) does not of course 
exclude the existence of more massive "excited" (and therefore in principle 
ultimately unstable) equilibrium states whose support depends on having 
more than minimal energy in thermal or other forms (such as that involved 
in differential rotation) which in the long run are subject to dissipation and 
loss by radiation but which in practice may be preserved over astrophysically 
or even cosmologically long time scales. The most important examples of 
such excited states are of course ordinary hot stars (including most notably 
those of the main, i.e. hydrogen burning, sequence) which are characteris- 
able by a mean central temperature say, in terms of which the pressure 
P will be given by the sum of a radiation contribution P ~ and of a 
non-relativistic gas contribution P ^ nQ where the relevant number density 
of non relativistic particles will be of the same order as the baryon number 
density n ^ p/mp. Evidently the radiation contribution will be dominant for 
Q'^ » n while the non relativistic gas contribution will be dominant for 
<< n. On substituting these formulae in the virial equilibrium condition 
(1.2) it can seen that the only criterion for radiation dominance is the mass 
M of the star, the critical value (which was first worked out by Eddington^^^l) 
being again given in order of magnitude by the inverse of the gravitational 
coupling constant m^^ i.e. by the same formula (6.1) as was obtained a few 
years later by Chandrasekhar for the more sharply definable upper mass limit 
for cold equilibrium : below this limit the dominant pressure contribution is 
that of the non-relativistic particles, whose substitution in the virial equi- 
librium condition (6.2) leads to a formula giving the characteristic density 
corresponding to a given characteristic central temprature in the form 




Fig. 1 - Logarithmic plot of characteristic orders of magnitude for density 
p against mass M relative to the standard values given by (1-1) and (1-6), 
which roughly characterise the Sun. 
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whereas when the stellar mass exceeds the critical value (1.6) it can be seen 
that it is the radiation gas pressure contribution that will be dominant so 
that instead of (1.9) one will obtain a less strongly mass dependent result 
expressible by 

The explanation for the fact that the typical masses of ordinary 
observable stars turn out to be comparable with the critical Eddington - 
Chadrasehkar mass value given by (1.6) is to be found in terms of the criterion 
for stability with respect to adiabatic variations, in which the pressure will 
vary as a function of density according to an approximately polytropic law of 
the form P oc where the index value is F = 5/3 in the non-relativistic limit 
to which (1.9) applies, but where F = 4/3 in the radiation dominated limit 
to which (1.10) applies. It is immediately obvious from the form of the virial 
condition (1.2) that it is necessary for stability that the effective polytropic 
index should exceed the critical value 4/3, i.e. precisely the same value that 
is characterises adiabatic perturbations of a radiation dominated gas. This 
means that stellar configurations above the critical mass (1.2) can at best 
be stabilised only marginally by their relatively small non-relativistic gas 
pressure contribution, and that for mass values a long way (more than two 
powers of ten) above the critical value M ^ l/m^^ stable equilibrium will in 
practice be extremely difficult to achieve. One would therefore expect that (in 
accordance with what is actually observed) formation of stars by gravitational 
condensation (with central heating according to the law 6 oc p^^^ that is 
obtained from both (1.9) and (1.10)) from initially diffuse gas clouds would 
inevitably produce objects below or not too far above the critical mass (1.6). 

As well as being limited in mass the conceivable range for ordi- 
nary stellar equilibrium states is of course also limited in temperature, which 
must exceed the Rydberg energy value, ~ e^rrie, that is the threshold for 
the ionisation of the gas that accounts for the opacity needed to delay the 
radiation loss of the thermal energy : on substitution of this minimal Ryd- 
berg temperature in (1.9) and comparison with the white dwarf equilibrium 
condition (1.4) it can be seen (see figure 1) that the smallest possible mass 
for an ionised stellar configuration with thermal pressure support is the same 
as the maximum possible value, M ^ e^/rrip^ for a cold planetary configura- 
tion, i.e. about the mass of Jupiter which is situated just at the lower end of 
the cold white dwarf range. At the opposite extreme the upper cut off to the 
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conceivable range of temperatures for ordinary stellar configurations is given 
by the electron positron pair creation temperature G ~ mg beyond which 
there is longer any possibility of stabilisation by a non relativistic electron 
gas contribution : on substitution of this pair creation temperature in (1.10) 
and comparison with the white dwarf equilibrium condition (1.4) it can be 
seen (see figure 1) that the highest possible characteristic density p for an 
ionised stellar configuration with thermal pressure support is obtained for a 
mass of the order of the Chandrasekar limit value (1.6) and is the same as 
the maximum possible value, p ~ me^mp, that is obtained at the upper end 
of the white dwarf range. 

Although energy loss by radiation from the outer ( "chromospheric" ) 
surface layers prevents them from lasting indefinitely, the "excited" stellar 
equilibrium states in the range delimited by the considerations of the preceed- 
ing paragraphs can nevertheless survive over astrophysically long timescales 
whose minimum value is determined by the minimal opacity contribution 
that results from Thompson scattering of photons by electrons with effec- 
tive cross section given in order of magnitude by cr ^ (e^/mg)^, which leads 
to a minimal evolution timescale r, for stars in the radiative mass range 
M >> l/rup'^ to which (1.10) applies, that will be given by 

r^^^, (1.11) 
TTif, rrip 

where e is the efficiency of conversion of rest mass into thermal energy by 
nuclear reactions. The most efficient thermonuclear energy production pro- 
cess is of course hydrogen burning which yeilds almost one per cent, e ~ 10^^ 
at a "main sequence" temperature at which stars spend most of the life- 
time allowed by (1.11), which works out at about 10^ years. For smaller stars 
with masses near or below the Eddington - Chandrasekhar critical value 
(1.6), other mechanisms come into play which increase the opacity and di- 
minish the rate of energy loss by radiation, giving timescales that for the 
smallest main sequence stars can greatly exceed even the present age of the 
universe which is of the order of 10^° years. The value of the relevant main 
sequence central temperature is derivable (by consideration of the probabil- 
ity of coulomb barrier tunnelling by the ionic reactants) as given in order 
of magnitude by the protonic analogue of the electronic Rydberg energy, i.e. 
^ e'^TTLp which is logarithmic between the minimal (ordinary electronic) 
Rydberg ionisation temperature ~ e^rrig. and of the maximal pair creation 
temperature ~ mg. 
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The foregoing considerations lead to the prediction (in full agree- 
ment with observation) not only that formation of stars in that mass range 
just above the the Eddington - Chandrasekhar critical value should have 
been relatively common, but also that most such moderately massive stars 
should already have passed the ends of the thermonuclear lifetimes and so 
been already obliged to face the issue of runaway gravitational collapse to 
densities in excess of the critical Michell Laplace limit '^^It^^l value 

(1.12) 

beyond which any description in Newtonian terms must be expected to break 
down, the usual formula for the scalar gravitational potential </? ^ M/R 
with R [M/pyf^ giving a result greater than unity, meaning that the 
gravitational energy is greater than the rest mass energy and hence that the 
escape velocity is greater than the speed of light, a situation that corresponds, 
in the General Relativistic formulation described in the following sections, 
to the light trapping mechanism that is the essence of the phenomenon that 
is commonly referred to as the formation of a black hole. Since the speed 
of light is normally supposed to represent an upper bound on the rate of 
propagation of causal influences of any kind, the infalling matter within the 
"horizon" (that is presumed to define the boundary of the region from which 
no light escapes) will become causally decoupled from the outside region, 
which thereby aquires the freedom to attain an equilibrium state of a new, 
essentially non - material "black hole" type, whose investigation will be the 
subject of the discussion in the following sections. Assuming that the density 
would retain its usual order of magnitude, i.e. that given by (1.1) Michell 
and Laplace estimated that light trapping would require a minimum mass 
of the order of 10^ times that of the sun (a value so gigantic that it was 
not taken seriously until, following up a suggestion by Lynden-Bell^^^l, its 
potential relevance to exotic quasar type phenomena in active galactic nuclei 
was pointed out by Hills'^^^ who noticed that it represents a threshold value 
for tidal disruption of ordinary starsl^^l). 

The paradox is that black holes in the relatively moderate mass 
range on which (by allowing for compressibility) our attention has so natu- 
rally been focussed by the line of astrophysical reasonning developped above 
i.e. a few times the Chandrasekhar limit (1.6), would appear in practice to be 
very rare, despite the high (predicted and observed) abundance of potential 
precursor stars. If these very common massive main sequence stars do not 
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become black holes, what happens instead ? 

One much discussed idea which emphatically does not give the 
correct explanation is that at the more than nuclear energy density that is 
attained at the upper limit, p k, rrip'^, of the neutron star density range, the 
equation of state deviates from the form (characterised by (1.7) and (1.5) ) on 
which the above reasonning is based in such a way as to allow the existence of 
ordinary material equilibrium states above the Chandrasekhar mass after all. 
The theoretical objection to this idea is that it would require the pressure to 
increase with density at a rate that would be incompatible with the causality 
condition that presumeably requires the corresponding sound (compression 
wave) speed with square given by dP/dp to be less than unity. Since P is 
small compared with p in the physically well understood low density regime, 
respect for this causality requirement that we should have P < p throughout 
the entire range so that satisfaction of the virial equilibrium condition for 
a given mass M entails that the density p cannot exceed the Michell limit 
value given by (1.12). This consideration does not rule ou the possibility 
of exotic ultra dense (e.g. quark nugget type) material equilibrium states 
with p » mp'^ and correspondingly with M << l/rrip^, but it does rule 
out their existence for higher mass values. A devil's advocate might still try 
to argue that one could still get round Chandrasekhar's upper mass limit 
by postulating some appropriately unorthodox relativistic gravitation theory 
for which the virial condition (1.2) itself would be suitably modified, but 
such theoretical gymnastics would seem to be pointless in view of the the 
complete absence of the slightest shred of observational evidence in favour of 
the existence of any such weird states as would be produced that way. The 
conclusion to be drawn from the intense astronomical activity of recent years 
is not just that plausible black hole candidates in the mass range just above 
the Chandrasekhar limit are comparitively rare, but also that there is no sign 
whatsoever of any alternative non-black hole type of cold (as opposed to hot 
stellar) equilibrium state at all in this mass range. 

We thus get back to the basic question of what actually has hap- 
penned to the numerous stars in the moderately massive range M > Xjm^ 
that due to of the comparitive shortness of the timescale (1.11) must have 
already burned out by now. The answer, which is implicit in the physics de- 
scribed in the preceeding paragraphs, can be presented in terms of several 
successive steps. To start with, since they are never far from instability, the 
radiation dominated stars inquestion will always tend to lose matter from 
their surface in the form of an outgoing winds which can carry away a very 
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significant fraction of he original mass during the last stages of the thermonu- 
clear lifetime. Secondly the dense burned out material that will accumulate 
in the core of the star will ultimately tend to evolve on its own almost in- 
dependently of the comparitively diffuse (even if much more massive) outer 
envelope layers. As soon as there is a degenerate central core in excess of the 
Chandrasekhar limit it can be expected to collapse by itself without waiting 
for the outer layers to be ready to follow. Surprised and shocked, these outer 
layers will thus be vulnerable to being blown away in a supernova type explo- 
sion by the energy released by the core collapse. The fact, refered to above, 
that the neutron star mass limit is rather larger, perhaps about double, that 
for the degenerate electron supported core, means that the core collapse can 
be expected to be halted, with formation of the shock that acts back on the 
outer layers, when the central density reaches that of neutron star matter. 
The conclusion (which is of course supported by a large amount of detailed 
numerical calculations by many workers) is that while small main sequence 
stars can obviously be expected to end up in white dwarf states, ones that 
are initially much more massive can be expected to end up by forming only 
slightly more massive neutron star remnants, the remainder of the mass be- 
ing dispersed in the form of a continous wind followed by an explosive burst, 
a picture whose broad outline is fully consistent with what is actually ob- 
served. Although the details are complicated and still highly controversial, 
it is easy to see that formation of a black hole is likely to be of more excep- 
tional occurrence, due to partial failure in a restricted parameter range of the 
supernova mass ejection process, or to subsequent accretion from a binary 
partner onto the neutron star remnant. 

2 The Example of Spherical Collapse. 

The collection of more or less well defined and physically plausible 
qualitative notions - unstoppable collapse with formation of an event horizon 
hiding the ensuing singularities - that constitute what may be referred to as 
the black hole paradigm was originally derived from the relatively tractable 
example provided by the spherically symmetric case, whose analysis will be 
the subject of the present section. In the case of actual equilibrium states, a 
considerable amount is now known about the non spherical generalisations 
that will be the main subject of later sections, but as far as dynamical evolu- 
tion is concerned, although we can draw a few general qualitative conclusions 
(such as the Hawking area theorem to be described in the next section) there 
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is still very little quantitative knowledge about what happens beyond the 
immediate neighbourghood of spherical symmetry. This makes it necessary 
to rely rather heavily on the spherical example despite of (and indeed even 
because of) the fact that even today it is still far from clear to what extent 
the lessons provided by the spherical model are generically valid. 

The mathematical analysis in this and all the following sections will 
be based on the standard Einstein theory of gravity as formulated in terms 
of a spacetime manifold with local coordinates x'^ say (/x = 0,1,2,3) and 
Lorentz signature metric field gfj_u (used for index lowering) that is governed 
by dynamical equations of the form 

G^" = R'"' - ^Rg^"" = SttT'"' (2.1) 

where Rfj_u (with trace R = R^^i) is the Ricci tensor of the spacetime metric 
g^^i, (see the appendix for definitions and notation conventions) and T^'^ is an 
appropriately chosen stress-momentum-energy density tensor whose form will 
depend on the kind of matter under consideration but which, for consistency 
with (2.1), must of course must always obey the "covariant conservation" law 

V^T^^ = (2.2) 

where is the standard operator of Riemannian covariant differentiation 
as defined with respect to g^y. 

As explained in the appendix we shall use an underline whenever nec- 
essary to distinguish quantities defined with respect to the geometry of an 
imbedded surface under consideration from the analogous quantities as de- 
fined with respect to the background geometry. As far as this present section 
is concerned the relevant imbedded surfaces are to be understood as con- 
sisting of the congruence of compact spacelike 2-surfaces generated by the 
spherical symmetry action, whose intrinsic Ricci curvature scalar will there- 
for, in accordance with this convention, be denoted by R to distinguish it 
from the background Ricci curvature scalar R. This allows the specification 
of what we shall refer to as the Misner Sharp mass function, M", by the 
formula 

{2Kf''^M^ = 2R- K^K'' (2.3) 

where K^^ is the extrinsic curvature vector of the spacelike two-surface, as 
defined in the appendix. This definition (whose right hand side is propor- 
tional to the mean of Christodoulou's mutually conjugate "mass aspect" 
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functions'^^l'^^l) has the advantage of being manifestly covariant and giving 
a result that is well defined as a strictly local field for arbitrary (not nec- 
essarily spherical) spacelike two surfaces, (in contrast with the related but 
only semi-local Hawking mass^^^', which involves surface integration over the 
two-surface, but which, like (2.3), was chosen so as to agree with the original 
mass specification of Misner and SharpP^l in the spherical limit with which 
we are concerned here). 

The specially convenient feature of the scalar field defined by (2.3) is 
that in the spherically symmetric case its derivative is directly related to the 
Einstein tensor of the gravitational field equations by (using square brackets 
to denote antisymmetrisation) the identity 

(2E')'/'V^M« = 2i^[,i^^] , h,, = i/i/Gp. , (2.4) 

using ± to indicate surface orthogonal projection (so that ^^j, = 9i-Lu—'g 
where g^^ is the fundamental tensor of the spherical two surfaces - see ap- 
pendix). This identity was first derived (though not in such a manifestly 
covariant form) by Misner and Sharp^^^l. It can be given a rather more ex- 
plicit form by introducing the usual circumferencial radius function r of the 
spheres, in terms of which their extrinsic curvature vector and Ricci scalar 
R (whose surface integral is of course Svr by the Gauss Bonnet identity) work 
out to be given simply by 

K, = -^V,r, i?=l. (2.5) 

so that (2.4) reduces to the form 

V^M» = r^G^^Vi^Vp] r, (2.6) 

whose derivation will now be described. 

One of the convenient features of a spherical (as opposed to more 
general) spacetime geometry is that it is possible to describe it in terms of 
an orthonormal tetrad of covectors 9^^, A = 0, 1, 2, 3, that is fully integrable 
in the sense that each one is proportional to the gradient of a corresponding 
preferred coordinate. These may be taken to be a provisionally unspecified 
space and time coordinate, and together with the usual spherical angle 
coordinates 9 and so that (using brackets to distinguish frame indices from 
coordinate indices) one has 
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e^^^dx^ = rd9, = r sin^ # , (2.7) 

of which the first pair of "outer" frame vectors are orthogonal to the two- 
spheres generated by the symmetry action while the last pair of "inner" 
frame vactors are tangential to them. The spherical symmetry is expressed 
by the condition that the three unknown metric coefficients r, ^p^, ^p-^ are 
all functions of x° and only. In terms of these quantities and of the fixed 
Minkowski frame metric (7^^ with signature (-1,1,1,1) the metric form 

ds"^ = gA<s>9\9\ dx^ dx" = -^^^dx^dx" + ^^^dx^dx" (2.8) 
where the "outer" part is given by 

-^i^Jx^'dx" = -<p^^dx" 2 + <p,'^dx^ 2 (2.9) 

while the metric within each two-sphere is given by the standard expression 

"^^Jx^dx" = r^{d9^ + sin^^ rf^ ^) . (2.10) 

Although the method most commonly given in textbooks proceeds 
by working out all the (forty) Christoffel components, the quickest way of 
evaluating the curvature tensor of a metric such as this is to use the Cartan 
techniquet^°l of proceeding via the calculation of the connection forms i^/i^^ 
which are got by solving the equations 

Vl,0\]=9-l,w,]\ , =0, (2.11) 

(using square and round brackets to indicate index antisymmetrisation and 
symmetrisation respectively, and with the understanding that the fixed Minkowski 
metric is used for lowering and raising of frame indices), the trick being that 
due to the antisymmetrisation there is no need to know the Christoffel con- 
nection components to carry out the covariant differentiation operations, the 
result being obtainable simply by replacing the covariant differentiation oper- 
ators by the corresponding partial differentiation operators 9^. The next 
step is to use this same (exterior differentiation) trick again in evaluating the 
corresponding curvature form 

Rt.u\ = + 2ti7j/eti7,]^^ . (2.12) 

It will be convenient for what follows to use a systematic shorthand no- 
tation whereby a suffix (q) or (i) is used to indicate the effect of differentiation 
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with respect to proper length in the space or time direction respectively, i.e. 
to indicate the corresponding frame components of the the covariant deriva- 
tive when acting on a scalar, so that in particular, for the circumferential 
radius function r itself we have 

1 dr 1 dr 

""t") " ' = —Q^ ■ (2.13) 

In terms of this notation scheme the six independent connection forms are 
found from (2.15) to be expressible as 

„ (0) _ MO) , am (2) _ _^^fl(3) 

= --fO''^ ^ ^.^^3) = --f(^'\ ■ (2-14) 

The Cartan formula (2.16) can now be used for the direct evaluation 
of the tetrad components Rq^^^ of the Rieman tensor, the only ones that 
are independent (bearing in mind that the spherical symmetry ensures that 
they are invariant under interchange of the indices (2) and (3)) being 

D (0) = ^Hom _ ^owm p (2) =^fi,^2_ 2 

{0){1) (1) ,„ ,„ ' -"'(2)(3) (3) »,2 V (0) (1) 

P (0) = ^(0)(0) _ ^(1) '^O(l) p (1) ^ ^(0) '^i(o) _ 

(0)(2) (2) ^ r If ' (1)(2) (2) ^ ^ ^ ' 

p (0) ^ ''^(o)(i) ^(1) '^i(o) _ ^(i)(o) ^(0) '^o(i) (9 1^] 

^(i)(2) (2) ^ r if, ~ r r if, ' ^ ' 

The corresponding frame components of the Einstein tensor will be given in 
terms of these by 

G =2R w 4- R (2) G = -2R (0) - R (2) 

'-^(0)(0) ^-"'(1)(2) (2) ^ -"'(2)(3) (3) ' ''(0)(2) (2) -"'(2){3) (3) ' 

'-^(0)(1) ^-"'(1){2) (2) ' ^(0){2) ^(1)(2) ^(2)(3) ^ ' 

G = -R (") - R - R (2 16) 

^{2){2) ''(0){1) (1) -"-(0X2) (2) -'^1)(2) (2) " \^^.±KJJ 

The resulting system can be considerably simplified by imposing that 
the coordinates be comoving with respect to the flow congruence determined 
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by the eigenvector of the energy momentum tensor, which is equivalent to 
the condition that the frame be such as to diagonahse the Einstein tensor, 
i.e. 

Cf(om =0 ^ '^I'^wd) = '^(d'^ko) ' (2-17) 
Subject to this requirement, which it is to be emphasized is not a physical re- 
striction but just a gauge condition, the only two Einstein tensor components 
still needed for the gravitational field equations can be seen to be expressible 
directly in terms of the Misner Sharp mass function specified by (2.3) whose 
explicit form is 

^" = ^(l + ^,o,'-^a,') ' (2-18) 

as 

_ W (o)_ (n TQ\ 

(0)(0) 2^ ' ^(i)(i) ^2^ ■ {^.J-JJ 

(1) (0) 

This can be seen to be just the frame component translation of the covariant 
version (2.6) of the Misner Sharp identity, whose derivation is thus completed. 
It is to be remarked that G^^^^^^^, the only remaining independent Einstein 
tensor component not given by(2.6), is not needed, because the equation in 
which it is involved will automatically hold as an identity whenever the other 
Einstein equations and the consistency condition (2.2) are satisfied. 

In considering the contributions to the right hand side of the Einstein 
equations (2.1) it is commonly convenient to work with a decomposition of 
the form 

T^^^ = T,/^ + T^'^'^ (2.20) 

in which T^^^'^^ is a "strictly material" contribution and Tp'^^ is an electro- 
magnetic field contribution given by 

T/" = — (r^'pF'^f - -Fp„Ff"'g>'A (2.21) 
47r V 4 / 

in terms of an electromagnetic gauge curvature field 

F^, = 2V[^A,] (2.22) 

(where square brackets indicate antisymmetrisation of the included indices) 
with (again necessarily) conserved source current. 



= -^VuF^" , V.J" = . (2.23) 

4:71 
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This formulation makes it possible to characterise the important "electrovac" 
case as that in which the source contributions T^j^^" and both vanish, the 
strict vacuum case being that in which the field F^j^ also vanishes. 

For many purposes, including those of the present section, it is suffi- 
cient to use a treatment in which the source contributions are not necessarily 
restricted to vanish but in which they are postulated to have the particularly 
simple form describable as that of a non conducting perfect fiuid, meaning 
that there is a preferred timelike unit vector and associated orthogonal 
projection tensor 'j^^, as characterised by 

If,,. = g^Lu + u^u^ , u^u" = -1 (2.24) 

with respect to which the material and electromagnetic source fields are spa- 
cially isotropic, meaning that they satisfy 

T^/'^ = pu^'u'' + P-f^"" , JK"^ = , (2.25) 

where the eigenvalues p and P are to be interpreted as the precise local values 
of the mass density and pressure whose characteristic mean values were the 
subject of discussion in the crude Newtonian order of magnitude treatment 
of the previous section. 

It will be sufficient for the purpose of the present section to further 
restrict our attention to the case of an adiabatic model, for which P is de- 
termined as a (not necessarily uniform) function of p allong each world line. 
As far as spherical applications are conserned there will be no further loss 
of generality in taking the model to be characterised by a pair of conserved 
number currents 

s>' = nu'' , n>' = nu'' , V^s'^ = , V^n'^ = (2.26) 

with n interpretable as the baryon number density and s as an entropy den- 
sity, in terms of which the mass density p is specified by a (uniform) equation 
of state function whose derivatives, interpretable as the effective temperature 
say and the effective mass-energy per particle or chemical potential p say, 
determine the corresponding pressure function P by the familiar relation 

P = se + np-p , e = ^, p=^ . (2.27) 

as an 

In terms of the corresponding thermal and particle four- momentum covectors, 

= , pp = pUp , (2.28) 
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and subject to the conservation laws (2.26) the perfect fluid equations of 
motion obtained from (2.2) are expressible just as the momentum transport 
equation 

2uP{sW^pQ,^ + nV[p^,]) = F^pJP . (2.29) 

From a computational point of view, this latter formulation'^^1'^^1'^^1 has the 
advantage (as compared with (2.2)) that, as in (2.11) and (2.12), the anti- 
symmetrised "exterior" nature of the derivation involved makes it possible to 
work it out directly, by direct substitution of the partial differentiation oper- 
ator = d/dx^ in place of the Riemannian operation V^, thereby making 
it possible to avoid having to go to the trouble of working out the Christoffel 
connection components. 

As a prerequisite to applying the formulae (2.19) in this perfect fluid 
case, it is necessary to impose the gauge restriction (2.17) to the effect that 
the coordinates should be comoving, which means that the timelike frame 
vector ^(0)^ is to be identified with the unit flow vector of the fluid as 
introduced in (2.24), so that by (2.25) and (2.26) the particle number and 
electric source currents will be given by 

JA* = en^" , = nO^"^ , (2.30) 

where the parameter e represents the electric charge per particle. It follows 
that the particle, entropy, and charge conservation laws (2.26) and (2.23) will 
be expressible simply by 

KVJ(0)=0, ©(0=0' e(o)=0 (2.31) 
or equivalently by 

iV^„, =0, = 0, Q(„, =0, (2.32) 

where A^, S and Q are fields respectively representing the total particle num- 
ber, entropy, and electric charge within the the corresponding sphere, which 
will be given as integrals over the interior of the sphere at a fixed value of 
the coordinate time x° by 

N f S f Q f 

Att J An J An J 

The corresponding frame components of the total energy momentum tensor 
are then then obtainable from from (2.21) and (2.25) as 

^(o)(o) = P + ^ ' ^(.)(.) = ^ - ^ ' (2.34) 
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where the appropriate electric field magnitude, E = F^_^^^^^, is obtainable by 
direct integration from the source equation (2.22) in the form 

E = ^ . (2.35) 

It can be seen from (2.21), by a similar integration, that the corresponding 
magnetic field component -^(2)(3) is necessarily zero, i.e. there can be no mag- 
netic monopole moment, on the assumption that (initially at least) there is 
a well behaved spherical centre from which the integrals in (2.33) are under- 
stood to be taken. It is to be remarked that the other (crossed) field com- 
ponents -^(o)(2)' -^(o)(3)' -^(i)(2)' -^(i)(3) vanish trivially as a locall requirement 
for spherical symmetry. 

Just as the classically familiar Coulombian form of the relation (2.35) 
is due to an judicious choice of definitions of the variables involved, so also the 
particularly astute Misner Sharp choice^^^' for the definition of the mass func- 
tion M" leads to a pseudo Newtonian form for the integral relation expressing 
the spacial constraint resulting from the first of the Einstein equations ob- 
tained from (2.19) which gives 



r E'^ 

= 471 r^{p + —)dr , (2.36) 
J Svr 

in which, as in (2.33) it is to be understood that the integral is taken over 
the interior of the relevant sphere at a fixed value of the comoving time 
coordinate, (starting from a central origin that is assumed to be regular, at 
least initially) which means that the radial variation will be expressible in 
terms of that of the space coordinate by the relation dr = r^^^^'^-^dx^ . 

Subject to the constraint (2.36) and the shift condition (2.17) which 
can be rewritten in the symmetrically equivalent form 

"Pormm = ^m^oii) ' (2-37) 
the only other Einstein equation that is needed is the dynamical equation 

^,0) =(^-4^^)^%) • (2-38) 

These two equations are to be solved in conjunction with the constraint 
obtained from the momentum transport equation (2.12 which takes the form 

{p + P)^^^^^=^,{Een-P^^^) (2.39) 
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in which the spacial pressure gradient is of course given by P^^^ = ri/i^jj +-50(1) • 
It is apparent at this stage that it will be convenient to introduce 
a modified, combined (electromagnetic as opposed to purely gravitational) 
mass function M say, given in terms of the original Misner Sharp mass func- 
tion M by 

n2 

M = + — (2.40) 
2r 

or equivalently by the more direct relation 

1 _ ^ = r^(K^K^ - E^) . (2.41) 
r 

It follows from the Misner Sharp identity (2.4) in conjunction with the Ein- 
stein equations (2.1) that the gradient of this combined mass function will 
be given in terms of the purely "material" contribution T^^'^'^ in the decom- 
position (2.20) by the manifestly covariant expression 

V^M = ^ V^g + 87rr2r,/,V[,V^]r . (2.42) 

For the case of a nonconducting perfect fiuid in the comoving frame, the 
constraint (2.36) and the dynamical equation (2.38) can be rewritten in terms 
of this new combined mass variable as 

M = Air f r^{p + neE—)dr , (2.43) 

^(1) 

and 

^(0) = -4vrPrV^„, . (2.44) 

of which the latter has the avantage of having the same simple form as that to 
which the original version (2.38) would reduce if no electromagnetic effects 
were present. An analogous remark applies also to (2.44) in any external 
region where the charge density ne vanishes. 

Whichever formulation is used, the solution of the system will in gen- 
eral require numerical computation, as in the pionneering attempt at an 
astrophysically realistic calculation by May and White^^^' or the important 
investigation of the possibility of naked singularity formation by Eardley and 
Smarr'^^]. It is however possible to obtain analytic solutions in special cases 
of which the most obvious are those in which the circumstances are such that 
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the right hand side of the constraint equation (2.39), and hence also that of 
(2.37) is zero so that we obtain 

V'o = 1 , ^,.)(o) = , (2.45) 

for a suitable normalisation of the time coordinate x° which in this particular 
case is adjustable to agree with the proper time allong the flow lines. Such 
a possibility obviously occurs in layers of matter that are uncharged and for 
which the pressure gradient P^^^ is zero, either because the configuration is 
homogeneous as in the classic prototype collapse calculation of Oppenheimer 
and Snyder'^^l or because the matter has a pressure free (so called "dust") 
equation of state. A fluid of this latter uncharged dust type (for which the 
flow will simply be geodesic) is characterisable by 

e = , p = mn , m^^^ = (2.46) 

where m is a constant mass per particle, so that in a layer of this type Q 
will be constant not just in time but also in space, while the combined mass 
function M will at least be constant in time : 



Q,,, = , M^o) = • (2-47) 

This means that the radial evolution equation , which by the definition of M 
will allways have the form 

9 9 2M g2 

will in this case be independently integrable for each flow line, since by (2.32), 
(2.45), and (2.47) the quantities Q, r^^^, and M appearing on the right will 
all be constants allong each separate flow line. The simplest possibility is 
the "parabolic" case corresponding to zero radial velocity in the large radius 
limit which is got by taking 

(Mr + Q^)J2Mr - Q'' 

r,^^ = 1 , x" - ^ — , c° = 2.49) 

(1) ' 3j\^2 ' (0) V ; 

where c° like M is an initially arbitrary constant along each flow line, i.e. a 
function only of . This latter comoving space variable can now be replaced 
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(except in the special case for which both c° and M are spacially uniform) 
by r in the "outer" part of the metric (2.8) which thereby acquires the form 



I 2M / 2M 02 

■^^Jx^dx'' = -(!- — + ^)(c/x«)2 + 2W(— - -^)dx''dr + dr^ (2.50) 



with the combined mass variable M now determined implicitly through its 
functional dependence on c° by the relation (2.49). 

The class of solutions specified by (2.49) and (2.50) is by no means 
simple and it is only comparitively recently (with the work of Eardley and 
Smarr'^^1 on the uncharged Q = case) that they have started to be exam- 
ined seriously from the point of view of questions such as naked singularity 
formation. They do however include the genuinely simple electrovac case for 
which the mass coefficient m in (2.46) is set equal to zero, which implies the 
constancy in space as well as comoving time of the combined mass variable 
M (but therefore not of the original Misner Sharp mass variable M" except 
in the Q = case for which they coincide). In this special electrovac case, as 
characterised by 

m = , M^^j = , (2.51) 

the flow just represents a geodesic test particle congruence, so there in fa- 
vor of the radius variable r so that there is no further loss of generality in 
imposing the parabolicity condition (2.49). This means that the form (2.50) 
with not just Q but now also M taken to be constant in space as well as time 
represents the most general spherical electrovac solution ( appart from the 
exceptional Robinson-Bertotti case for which c° is uniform, so that'^^^ one 
obtains a tubular universe with constant radius r throughout). This solution 
can be seen to be automatically stationary since all dependence on x° has 
dropped out, which in the spherical case means more particularly that it 
must be static, i.e. that it is unaffected not only by displacements but also 
by reversals of a certain preferred time coordinate, t say, that is determined 
(modulo a constant of integration) by the differential relation 



dt = dx'^ - , ^ . (2.52) 



V2Mr - 
r2 - 2Mr + 

Replacement of the comoving proper time coordinate x° by this preferred 
time coordinate t leads to the manifestly static form 




^.Jx^dx"" = -{l- — + ^\de — + dr' (2.53) 
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that was originally derived by Reissner and Nordstrom on the basis of the 
postulate of staticity at the outset. The present approach, showing how static- 
ity is obtained automatically as a consequence of spherical symmetry in the 
source free case, amounts to a demonstration of what is known as Birkhoff 's 
theorem. 

Our parabollically infalling version (2.50) has the significant advan- 
tage over the algebraically simpler historic form (2.53) that it remains well 
behaved on the "Killing horizons" '^^1 , i.e. stationary null hypersurfaces that 
occur, whenever < M^, at the roots 

r = r± , r^ = -M ± ^ AP - (2.54) 

whereas the manifestly static version (2.53) is singular there. However al- 
though it is sufficient for describing the outside of a collapsing spherical 
charged or neutral star model, even the more sophisticated version (2.50) 
has the limitation of being geodesically incomplete even when extended over 
the full coordinate range < r < oo, — oo < x° < oo. The geometrically 
complete manifold was first described in the pure vacuum case, Q = O by 
Kruskal and Szekeres^^^^ and in the generic case < by Graves and 
gj^jjj[38]^ while for the special "maximally charged" limit case = the 
corresponding construction was first carried out rather later by myself'^^' . 

It was for the purpose of describing such extensions that I first intro- 
duced the representational technique of conformal projection (the space 
time analogue of the beloved Mercator projection of terrestrial navigators) 
that has since been generally adopted as a standard tool for understanding 
the topological and causal structure of any timelike two dimensional manifold 
or submanifold, the idea being to first convert the metric into null coordinate 
form 

~g^^dx''dx'' = -'^du+du- (2.55) 

which is always locally possible for some conformal factor determined as 
a function of the null coordinates u+ and u^, and then to take advantage of 
the fact that this null form is preserved, only the functional dependence of 
the conformal factor being altered, by a transformation u+ i— > n+, n_ i-^ £t-, 
v[/ hH. \1/ whereby each of the null cordinates is replaced by an arbitrary 
function only of itself, which one is free to choose in such a way as to cover 
what from a metric or affine point of view might be an infinite region by a 
finite coordinate range which can thus be plotted directly as a diagram. 



25 



Fig. 2 - Facsimiles of the original CP. diagrams^^^^^'^'^^ for the limiting cases 
= 0, i.e. Kruskal-Schwarzschild (including disconnected branch represent- 
ing analyutic extension to region r < {)), and = M"^ , i.e. "maximally 
charged" Reissner- Nordstrom. 
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The use of this method of representation of event horizons in two di- 
mensional manifolds was suggested to me by the example from a related but 
rather different context (namely the study of the distant assymptotically flat 
outer regions) of Penrose's conformal boundary procedure procedure'^^1 (as 
used for setting up what is known to the initiated by the term prounounced 
as "Scri"). However the concepts of a conformal boundary and a conformal 
projection should not be confused. The Penrose conformal boundary (Scri) 
concept is not restricted to two dimensions, but on the other hand it is de- 
pendent on rather severe assymptotic regularity conditions that may fail in 
many relevant applications, whereas the conformal projection (CP.) tech- 
nique is limited to two dimensional sections but not otherwise restricted, so 
that in particular it is very useful for the analysis of boundaries that may be 
singular. 

In the present case the most obvious choice to start with is to take the 
the incoming and outgoing stationarity preserving null coordinates which are 
defined (modulo an arbitrary choice of origin ) by 

t'^ diT 

du^ =dt± — — (2.56) 

r2 - 2Mr + ^ ' 

whose separate substitution gives the forms 

-g^Jxf'dx" = ^ {du^y ± du^dr . (2.57) 

These null inflowing and outflowing coordinate forms (like our original parabol- 
lically inflowing form (2.53)) are locally well behaved on the Killing horizons 
at r = r-t but nevertheless still incomplete. When both are substituted to- 
gether, with the radius variable now considered no longer as a coordinate in 
its own right but just as a function of the null coordinates, the metric aquires 
the required doubly null form (2.55) with the conformal factor given by 

r2 - 2Mr + 

* = -2 — (2-58) . 

In the generic case 

g2 < M^, the integration of (2.56) gives the explicit 

expressions 

= t ±r ± Inlr — r 1 I ± Inlr — r_| , (2.59) 

2k+ 2/t_ 
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from which, in terms of the decay constants of the Kilhng horizons which are 
given by 

K± = ± 5 , (2.60) 

r_i- 

one obtains the functional dependence of the variable r in (2.59) in the im- 
plicit form 

2r H ln|r — r+| ln|r — r_| = m+ — m" , (2.61) 

whose unambiguous solution requires the specification that r should lie in 
some particular one of the three possible ranges characterised by the condi- 
tion that neither, just one, or both of the quantities r—r± be positive. Except 
for the first of these three possibilities, which includes the value r = that 
correponds to an irremovable geometric singularity, the resulting conformal 
factor \l/ will be regular over the full coordinate range oo < u± < oo, but 
from the point of view of completeness the metric version given by (2.58) 
and (2.61) is no improvement on the traditional manifestly stationary ver- 
sion (2.53) : all that has been achieved is to push the Killing horizons out of 
the coordinate chart, but not to regularise them. However a genuine regular- 
isation is now easily obtainable by a conformal coordinate transformation to 
a a new null coordinate form 

-^^^dx^dx"" = ^du+du- . (2.62) 

Depending on whether it is an "outer" Killing horizon at r = r_|_ or an "inner" 
one at r = r_ that one wishes to cover, it suffices to take 

f = Yr''"7i . (2-63) 
where r is now given implicitly by 

±(r± - r)e^''^''\r - r^|'^±/'^T = u+u-. (2.64) 

This relation shows in particular how, at the outer horizon the decay param- 
eter is interpretable as measuring the exponential relation between the 
affine time parameter u- and the group parameter u- . 

The unambiguous solution of (2.64) requires the specification that r 
should lie in one or other of just two possible ranges characterised respectively 



28 



by r < Tzp and r > r±, and except for the irremovable geometric singularity 
at r = in the first of these two ranges, the new conformal factor ^ will 
be regular over the full coordinate range oo < u± < cxo, of the new null 
coordinates, including the locus r = r± which can be seen to consist of two 
intersecting Killing horizons characterised in the new coordinates by Ct+ = 
and u- = respectively. The stationarity group transported coordinates 
u+ and w of the original system in the more restricted patches on either 
side of these now regularized horizons are given in terms of the new ones, 
which can be seen to be characterised by the property of measuring affine 
distance allong the regularised horizon at r = r±, by relations of the simple 
exponential form 

= ^1-11+ 1 , K±u' = — Inlw"! . (2.64) 

The use of transformations of the simple form (2.64) allows us, ac- 
cording to choice, to cover either the locus r = r+ or the locus r = r„ with 
a regular coordinate chart, but it does not allow us to cover both at once. 
Nevertheless since the alternative kinds of chart overlapp (in the intermedi- 
ate range r_ < r < r+ where both are perfectly regular) they can be used as 
successive patches to build up a maximally estended manifold in the man- 
ner first described (in terms of a somewhat different system) by Graves and 
Brill [38]. It is for the purpose of visualising the final result of such successive 
extensions that the CP. (conformal projection) technique^'^^lt^^] is particu- 
larly useful. If one is willing to sacrifice the desideratum of having a simple 
analytic expression such as (2.64), there is no obstacle in principle to the 
introduction of further modified null coordinates u^saj whose range covers 
the entire maximally extended manifold. 

An invaluable practical feature of a CP. diagram of the kind obtained 
by plotting such coordinates directly (traditionally with a diagonal orienta- 
tion) on a flat screen or page (as in figures 2 and 3 ) is that, as far as the 
essential causal and topological features are concerned, it does not matter 
whether or not one knows the precise functional form of the functional re- 
lation between the original (restricted) and new (extended) null coordinate 
systems : provided the linear (diagonal) representation of the null congru- 
ences is preserved, any smooth (not necessarily analytic) deformations are 
admissible. This means that (provided it is not restricted by the inclusion of 
too much detail) any CP. diagram that has been constructed as a rough free- 
hand sketch has the beautiful feature of being interpretable post facto as an 
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Fig. 3 - Facsimile of the original CP. diagravr^'^^^ for the Graves - Brill ex- 
tension of Reissner Nordstrom for < < together with previously un- 
published "new look" version in which complete compactification is achieved 
by letting the scale for successive universes tend to zero at the extremities of 
the chain. 
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accurate representation in terms of null coordinates whose precise specifica- 
tion (if one were interested) could in principle by found out later by carrying 
out empirical measurements on the sketch. 
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3. Qualitative theory of non-spherical Black Hole formation. 

Whereas a considerable amount is known about non spherical 
black hole equilibrium states (to which the subsequent sections will be de- 
voted) as also about non stationary states of spherical collapse (the subject 
of the previous section) the subject of generic nonspherical gravitational col- 
lapse and black hole formation still consists mainly of a few vague, quali- 
tative, and for the most part far from rigourously established notions, that 
are largely inspired by the spherical example. The question of the extent 
to which various features of spherical collapse scenarios may be taken over 
to more general situations has long been and still remains a subject of ani- 
mated debate. The unreliability of the spherical example as a guide to more 
general cases is shown by the case of Birkhoff 's theorem, to the effect that 
(as was demonstrated in the previous section) the source free (strict or elec- 
tromagnetic) vacuum outside a collapsing spherical object must necessarily 
be static (i.e. not only time independent but even time reversal invariant) 
whereas in the non spherical case it need not even be stationary (i.e. time 
independent) in view of the possibility of gravitational and electromagnetic 
radiation whose absence, exceptionally, in the spherical case is due to the 
absence of any scalar part of either the electromagnetic field which is purely 
vectorial or the gravitational field which is purely (i.e. tracelessly) tensorial, 
at least in Einstein's theory to which our discussion here is restricted. 

Among the features that are generally thought to survive the the 
breaking of spherical symmetry of the collapse, some of the most important 
may be listed as follows : 

(1) The ultimate formation of a singularity of some kind ( not necesarily 
just a simple density singularity, but something sufficient to prevent affine 
completeness) was shown by the work of Penrose and Hawking to occur very 
generally , but its generic nature is still not well established. 

(2) The phenomenon for which Penrose'^^^ coined the term cosmic censor- 
ship, whereby the singularities are hidden from the outside asymptotically 
fiat universe behind a regular event horizon bounding the region for which 
Wheeler coined the term Black Hole would appear to be stable against mod- 
erate perturbations from spherical symetry and from the uniformity of the 
homogeneous Oppenheimer Snyder collapse scenario that provides its sim- 
plest example. Nevertheless much recent workt^^!'^''''^^'!^^!!'^^! has made it clear 
that sufficiently (one might be tempted to say unnaturally) large deviations 
from uniformity can bring about the occurrence of non-trivial naked singular- 
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Fig. 4 - Illustration of the standard black hole paradigm by a CP. diagram 
representing a two dimensional radial section through a collapsing body with 
cross hatched shading of the D.O.C. (domain of outer communications). 

ities, i.e. ones from which hght can escape to large asymptotic distances, so 
although the regular black hole scenario, as governed by the cosmic censorship 
postulate, may plausibly provide a generic description of astrophysically re- 
alistic collapses, its mathematical generality would seem to be more severely 
circumscribed than was once thought. 

(3) Although the vacuum region outside a generic collapsing body will not 
become immediately static (as it must, by Birkhoff 's theorem in the spherical 
case) it is nevertheless to be expected that the energy of non stationary 
oscillatons will ultimately be radiated away so that in the end the vaccuum 
region outside an (isolated) collapsing body will settle down asymptotically 
towards an ultimate equilibrium that is stationary at least in the weak sense 
of being invariant under the action of a Killing vector field that is timelike 
at large distances even if not everywhere outside the horizon. 

Experience with the Schwarzshild and Reissner Nordstrom exam- 
ples (as described in the previous section) shows however that whereas the 
physical collapse situation may be regular in the past, starting with an or- 
dinary well behaved assymptotically flat Cauchy initial value hypersurface, 
the asymptotically approached equilibrium metric may have a "white hole" 
region including singularities in the past, so the strongest regularity condi- 
tion it can be expected to satisfy is assymptotic predictability, meaning that 
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'caustic 



Fig. 5 - Illustration of possibility of cusp on past event horizon. It is shown 
how a null generator is attainable as limit of timelike lines escaping from 
sequence of points approaching horizon from outside. 



there exists a partial Cauchy surface (a not necessarily complete globally 
spacelike, i.e. achronal, hypersurface) extending in from outer infinity to the 
black hole horizon and at governing (i.e. intercepting all sufficiently extended 
past directed timelike lines from) not necessarily all of its future (as would be 
required for a strict Cauchy surface) but at least a part consisting of a regular 
asymptotically fiat domain of outer communications with inner bound on a 
regular black hole horizon. 

The future boundary of the region governed by the partial Cauchy 
surface (which in the standard Reissner Nordstrom example occurs at r = r_) 
is called its Cauchy horizon, and is an example of what is commonly described 
as a local "future event horizon", whereas a black hole horizon (which in 
the Reissner Nordstrom example occurs at r = r+), i.e. the boundary of the 
region from which a future directed timelike line can be extended to the outer 
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asymptotically flat region, is analogously describable as being a local "past 
event horizon". Both kinds obviously belong to the category of "achronal 
boundaries" , meaning boundaries of which no two points are connectable by a 
strictly timelike curves. Local future and past event horizons are characterised 
more particulary by the property (which was first systematically exploited 
by Penrose'^^l'^^l and Hawkingt^'^'lt^l) of being generated locally by by null 
geodesies with no respectively past or future end points (which means that 
they are ordinary null hypersurfaces wherever they are smooth, but that 
there null generators may reach caustics when extrapolated to respectively 
the future of the past). In the case of a black hole horizon this property (see 
figure) is not just local but global i.e. its null generators can never reach a 
future end pont no matter how far they are extrapolated. 

To draw quantitative conclusions from these considerations it is 
necessary to recapitulate some of the standard kinematic properties of gen- 
erating congruences. To start with we recall that for any vector field £^ has 
an (unnormalised) acceleration £^ that is related to its Lie derivative with 
respect to itself by 

ia^ = + 1 v^(r4) , t = rv.F . (3.1) 

In order for the field to be normal to a hypersurface, it must satisfy the 
Frobenius integrability condition 

^[^yjp] = , (3.2) 

which implies 

£[/££p] = r4V[/p] . (3.3) 

It is apparent from (3.3) that in the particular case of a hyper- 
surface that is locally null, i.e. whose normal satisfies 

r4 = , 

this normal must automatically satisfy the geodesic equation 

A/p] = . (3.4) 

Let us introduce a second null vector say, transverse to the null hypersur- 
face and normalised so that 

r4 = , r4 = -1 • (3.5) 
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Such a vector can be uniquely specified by the condition that it be normal 
to some given spacelike 2-surface S say in the horizon at the point under 
consideration, in which case the corresponding rank 2 projection tensor 

Tu = 9^u + i^t + i^L , ru = 2 (3.6) 

will be interpretable (see appendix) as the (first) fundamental tensor of the 
spacelike 2-surface. 

With respect to any previously specified normalisation, corre- 
sponding to a time parametrisaton such that = dx'^/dt, we can define 
a corresponding affine time parametrisation r say whose relation to the orig- 
inal time parameter r specifies a corresponding decay coefficient k, in terms 
of which the non affine geodesic equation (3.4) takes the form 

4 = (3-7) 

where explicitly 

K = (In f)- = f/i = -ti^ . (3.8) 

The usual way of defining the divergence 6 and the (automatically 
real) magnitude a of the (automatically symmetric) shear rate tensor a^;, of 
the null generators is via the projection 

Of.. = Wf.Wv^/a , e = e%, a^, = e^, - '-ff§^^ , = 2a^,a^'' . (3.9) 

Taking the contraction of the defining identity of the Riemann tensor, which 
for any field ^^ whatsoever gives the consequent identity 

^'^v^(v.r) = v,.(rv.r) - {vnv)vj^ - R^j^r, (3.10) 

one obtains, for the case of a congruence of null geodesic generators as char- 
acterised above, the famous Penrose t^^l null version 

e - kO = -'-{6^ + a^) - R^,Fr , (3.11) 

of the equation whose analogue for a timelike congruence was first brought 
to attention by Raychaudhuri'^^l , where a dot denotes differentiation with 
respect to an arbitrary time parametrisation, whose adjustment to be affine 
can be used to get rid of the k, term on the left hand side. The special impor- 
tance of this equation comes from the fact that one is then left with a right 
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hand side that is strictly non positive provided the Ricci tensor is determined 
by the Einstein equations (2.1) either for a vacuum or at least with an energy 
tensor T^'^ that, as is the case for all the usual macroscopic matter models, 
is such as to satisfy the appropriate energy inequality T^'^ifj^ii, = 0. 
The Penrose inequality 

r')- > ^ (3.12) 

that is obtained for the rate of affine variation under such conditions can be 
immediately used to see that if 6 is ever negative then there will inevitably 
be a caustic where it diverges to infinity within an affine distance bounded 
above by 2/|^| in the future. Penrose's original exploitation of this result was 
for the purpose of demonstrating the inevitability of some kind of singularity 
formation to thr future of any closed trapped surface on a well behaved initial 
hypersurface by showing the affine bounedness (which would be impossible 
in the absence of a singularity) of the future event horizon bounding the 
future of the closed trapped surface, where this term is understood to mean a 
compact spacelike topologically spherical 2-surface S for which the divergence 
6 of the null normals is everywhere negative. 

In terms of the background tensor curvature formalismt^^l de- 
scribed in the appendix, it can be seen (from (A9)) that the specifications 
(3.9) for the divergence and shear of the outgoing null congruence £^ can be 
rewritten in terms of the second fundamental tensor K^y'' (which is equiva- 
lent to what is referred to by Hartley and Tucker^^^l as the shape tensor) of 
the spacelike 2-surface S and of the corresponding curvature vector and 
trace free (and conformally invariant) conformation tensor C^,y^ as 

^ixu = —K^i^^ip , 9 = —K iu , (Jfiu = —Cfj^u'^ip . (3.13) 

Similarly for the ingoing null congruence as specified by £^ (for which we 
are still assuming the normalisation condition (3.5)) the corresponding di- 
vergence and shear will be given by 

0^.u = -K,/ip , e = -K^l, , ~a,, = -C,/lp . (3.14) 

The usual situation for an approximately spherical 2-surface at 
approximately constant time in an approximately fiat background is to have 
ingoing null normals that converge, ^ < but outgoing ones that diverge. 
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9 > so that the product, which can be seen from (3.6) to be expressible in 
manifestly the normalisation independent form 

99 = -\k,K' , (3.15) 

will be negative, whereas if 9 changes sign and becomes negative, with 9 still 
also negative, then the product will change sign also, i.e. the curvature vector 
will change from being spacelike to being timelike. A special interest 
applies to the marginally trapped case characterised by 6* = everywhere, for 
which Hawkingt"*^! has introduced the term apparent horizon. Evidently 
such a marginally trapped surface may be described as one for which the 
curvature vector A'^ is null. 

Following the Penrose application of (3.12) to the future of a 
closed trapped surface, Hawking[''^ll^°l[^^l pointed out that a very powerful 
result can be obtained by applying it to the black hole horizon itself, using 
the condition that the future generators of a black hole horizon can never 
terminate, which implies that the generators of such a horizon can never 
have negative divergence 9. Noting that 9 is interpretable as specifying the 
rate of variation of the measure of a 2-surface element dragged along by the 
generators according to the formula 

{dSy = 9dS , (3.16) 

and applying this to the integrated area 

A = jdS (3.17) 

of a 2 dimensional spacelike section through the black hole horizon. Hawking 
obtained the (now famous) law to the effect that the horizon must evolve 
with time according to the inequalities 

A> (j^9dS>0 , (3.18) 

the extra inequality on the right being to allow for the fact that in addition to 
the area increase resulting from smooth expansion there is also the possibility 
of an additional increase due to the branching off of new generators from a 
caustic (see diagram). In particular if two black holes with areas A i and A 2 
merge to form a combined black hole with area A 3 then we must have the 
strict inequality 

^3 >^ 1+^2- (3.19) 
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Before closing this section and going on to restrict our attention to 
states of stationary equilibrium, it is to be remarked that though the locality 
of a marginally trapped surface or "aparent horizon" , as characterised locally 
by K^K^ = 0, is merely an inner bound on the location of (teleologically 
defined) true horizon, it may nevertheless give a very good approximation 
to the localisation of the true horizon in the almost stationary limit when 
the situation is not too strongly different from its ultimate equilibrium state, 
in which case the approximate stationarity will determine a corresponding 
approximately well defined and in general non affine time parametrisation 
on the horizon, so that there will be a correspondingly well defined decay 
constant k = f/f where r is the corresponding proper time. Under such 
conditions the ratio between the values dS^ and dSi of the measure of a 
generator transported surface element dS between times to and ti can be 
seent^^lt^^l to be given approximately by an expression of the non-teleological 
form 

= r edt^—Ddt (3.20) 

OOq J Jto AC 

which gives the Hartle Hawking formula^^^' 

A ^ (f—dS , (3.21) 




K 

where the rate of effective dissipation is given by 

/^ = ^^^ + T'^%4. (3.22) 

The discovery'^ ''I'^^l of the laws (3.18) and (3.22) suggested an 
obvious thermodynamic analogy, with A proportional to the entropy and k to 
the temperature. The deeper significance of this analogy was first guessed by 
Beckenstein'^^1 and later established by the discovery of Hawking radiation'^'^1 , 
which is associated with a temperature given exactly by fi;/27r in Plank units, 
corresponding to an entropy given exactly by A/ 4:. The crux of the analogy is 
constituted by the "zeroth law" that will be established in the next section. 



39 



4. Rotating Equilibrium States : zamos and local properties of a 
Killing horizon. 

The concept of stationarity that is relevant to the theory of black hole 
equilibrium states means the condition the the spacetime is invariant under 
the action generated by an assymptotically timelike Killing vector 

|, ^i,K) = (4.1) 

It is to be noticed that this definition is slightly weaker than the one com- 
monly used in other contexts where it is stipulated tht the stationarity Killing 
vector be timelike not just at large asymptotic distances but throughout, 
which would exclude the existence of the "ergorgions" which are of impor- 
tance not just in black hole thory but even in the theory of ultrarapidly 
rotating stars. 

All that follows will be based on the postulate that 1 decided to adopt 
when 1 first looked into this area of work in the 1960's, namely that the 
stationary spacetime under consideration is also characterised by axisymetry, 
meaning that it is also invariant under the action generated by a spacelike 
Killing vector 

d 

m'^ ^ — , V(^m,) = (4.2) 

whose action is periodic, with closed circular (or, on the axis itself, fixed 
point) trajectories allong which the group parameter </> is therfore also pe- 
riodic, with period 2tt for the standard normalisation. Relaxation of the 
requirement that m'^ be spacelike is mathematicUy conceivable but physi- 
cally inappropriate since, in view of the periodicity, it obviously violate the 
causality requirement that there exist no closed timelike or null lines. In an 
asymptotically fiat background it is inevitablet^^^ that this second symmetry 
commutes with the first, i.e. 

rV^m'^ - m'^V^A;^ = 0. (4.3) 

Just as it is plausible that a isolated system with or without a cen- 
tral black hole should tend towards a stationary equilibrium state so also 
it is plausible, particularly in a context where gravitational radiation needs 
to be taken into account, that under natural conditions the stationary state 
would also have to be axisymmetric. It is of course possible mathematically 
to construct artificial counterexamples (such as the Dedekind ellipsoids to 
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which Chandrasekharf^^l has drawn attention), but under natural astrophys- 
ical conditions it is very hard to imagine stationary black hole scenarios for 
which the axisymmetry assumption would be in danger of failing. (For the 
case of of an isolated black hole with a vaccuum or simple gaseous exterior 
considerable effort has been invested, most notably by Hawkingt^^H^^lt^l, in 
attempts to prove that axisymmetry of equilibrium states is not just phys- 
ically plausible but mathematically necessary. However the crucial result, 
describable as Hawking's "Strong Rigidity Theorem", to this effect is itself 
based on a postulate of analyticity that is also physically plausible but whose 
mathematically justification requires assuming the axisymmetry one wanted 
to prove in the first place, so that as a mathematical challenge the problem 
remains wide open.) 

In any study of stationary axisymmetric systems an important role is 
played by the Killing vector invariants 

V = -kf'k^ , W = k^m^ , X = m^m^, (4.4) 

and by the determinant 

w'^ = VX\W'^ = ^JJ^yw""^ , w^^ = 2fc[^m^] . (4.5) 
and the ratio 

W 

w = -— , X > (4.6) 

which is well defined wherever X is strictly positive, which by the causality 
condition that m'^ is spacelike, will hold everwhere except on the symmetry 
axis itself where m'^ reduces to a zero vector so that X and W both vanish, 
making u undefined. 

Other contractions of interest are the energy, E say, and angular mo- 
mentum L say, of a particle with momentum covector p^, as given by 

E = -k'p, , L = Tjfp, , (4.7) 

which are of course conserved for free orbits : 

u''VuP^, = , = mu^ , u^'u^, = -1 ^ u^V^E = , u'^V^L = . 

(4.8) 

The quantities zu and u defined by (4.5) and (4.6) are of particular 
interest in the context of circular flow , meaning flow allong trajectories on 
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circles generated by the Killing vectors, i.e. with unit flow vector u'^ of the 
form 

= a{k^ + nm") (4.9) 

where the coefficient Q = d(f)/dt is the angular velocity of the trajectory. The 
important "Kepplerian" special class of circular trajectories consists of those 
that are free in the sense of (4.8), a possibility which typically will exist only 
in a restricted equatorial plane. A more generally defined class that is of 
more immediate (though mathematical rather than physical) interest for our 
present purposes consists of what Bardeent^^l has called "zamo" trajectories 
(short for zero angular momentum orbiters) which are characterised by the 
(obviously non Kepplerian) condition of having L = 0, which can be seen to 
be equivalent to the condition that their angular velocity be given directly 
hj Q = uj. The obvious interest of w as defined by (4.5) in this context is 
that its reality, i.e. the positivity condition 

>0 (4.10) 

is evidently the necessary and sufficient condition for the existence of a 
strictly timelike zamo at the position in question. The importance of this 
purely local condition is that subject to very weak hypotheses it can also 
be shown as described later on below, to characterize the domain 
of outer communications, whose (globally defined) boundary at the surfaces 
of the black hole region, will be characterisable (locally, this is what is so 
convenient) as a "zamosphere" , where w = 0, on which the zamo's become 
null. 

A more frequently discussed t^*^' but in the final instance less important 
analogue of the zamosphere is the "ergosphere" , where ^ = 0, i.e. on which 
the stationarity generator fc^ becomes null. The interest of this is that it 
bounds the "ergoregion" characterised hj V < within which a free parti- 
cle energy E, as defined by (4.7) can become negative, whereas outside the 
ergo region, i.e. wherever k'^ is timelike, the free particle energy is bounded 
below by the condition E > rrw/V . The existence of an ergoregion (un- 
less confined within the horizon as in non rotating case) makes possible the 
extraction of energy from the background by the mechanisms such as the 
Penrose process'^^' whereby a particle coming in with energy E\ splits into 
a part with negative energy, i?2 < 0, and an outgoing part with energy Ej, 
which, by conservation of the sum, E2 + E^ = Ei must exceed the initial 
energy, E^ > Ei. 
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For charged orbits, as given by 

u'^W^p^ = F^^u" , (4.11) 

in a stationary field, 

kCA^ = ^ F^^k" + = , $ = -k^A^ (4.12) 

one gets conservation not of tlie ordinary energy E but of a generalised 
gauge dependent generalisation £ say, constructed from the gauge dependent 
generalised momentum covector 

7r^=p^ + eA^, (4.13) 

i.e. one gets 

u''\/^£ = 0, £ = -h''TT^ = E + e^ , (4.14) 

This generalised evergy can be negative even outside of the ergosphere, where 
> so that the lower bound will be given by £^ > m\/V + e$. It is apparent 
that there will be an extended electric ergoregiont^^lt^^^t^^l characterised for a 
given charge to mass ratio by the possibhty of negative energy for at least 
one sign of the charge ±e, for which the condition is just 

V < (-)2$2 . (4.15) 
m 

The significance of this relation is of course dependent on how the energy 
is calibrated, the usual assymptotic specification being not the only one of 
interest : another possibility of particular interest '^^^t^^^t^^^ for the specialised 
theory of non rotating black holesI'^^t^^lt^^^t^^lt^'^Jt^^lt^^] is to calibrate with re- 
spect to the horizon which is possible in that case because of its uniform 
potential condition which will be demonstrated below. 

The properties of a stationary axisymmetric system simplify enor- 
mously under conditions of what I call circularity which in practice are al- 
most sure to be satisfied in the applications that are relevant to black holes. 
In the case of an electromagnetic source current j'' this means just that it 
should be a linear combination of the Killing vectors which is equivalent to 
the requirement 

jlA'tu^-l = , (4.16) 

while for a gravitational source it means inthe case of a perfect fluid just that 
the corresponding flow vector should satisfy the ananlogue of (16) which 
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is equivalent to the condition (4.9) for the corresponding flow trajectories to 
be simple circular orbits. For a more general material source (for which a 
preferred reference vector might not be deflned) circularity is to be under- 
stood as meaning that the relevant material energy and angular momentum 
flux vectors k^T^^ ^ and myT^" ^ should have the same property, this general 
circularity condition being expressible as 

kuTJ'^^Wp^^ = , myT^''[^Wp^] = . (4.17) 

This condition can in principle fail in a star that is partly solid'^^l (as 
in the case of a neutron star crust ) or even in a stricty perfect inviscid fluid 
star where there is convection, but for the more plausibly relevant case of 
viscous fluid the possibility of other than circular motion can be ruled out 
because it would inevitably produce thermal dissiption and thus violate the 
requirement of strict stationarity. 

The crucial simpliflcation that one gets in such circumstances is pro- 
vided by the circularity i/ieorem'^^l'^^l'^^l which (generalising a result flrst 
demonstrated in the case of a pure vaccum by Papapetrou'^*^' and for an un- 
charged perfect fluid by Kundt and Trumper'^^]) tells us that the system will 
be orthogonally transitive, meaning that the circular trajectories generated 
by the two Killing vectors will be orthogonal to a congruence of two dimensin- 
ional surfaces which (must obviously be spacelike where, and only where, the 
zamos are timelike) throughout any continuous region connected to the rota- 
tion axis within which the source circularity conditions (4.16) and (4.17) are 
satisfled. Thus if the source circularity conditions are satisfied everywhere 
one gets orthogonal transitivity everywhere. The well known Frobenius con- 
dition for such orthogonal transitivity is expressible as the requirement that 
the twist vectors 

= ie'^'^P-kyVpk^ , V = le^^'^myVprn, , (4.18) 
be orthogonal to the Killng vectors, i.e. that we should have 

uj^'m^ = , i^^'kf, = . (4.19) 

Since these contractions must both vanish identically on the axis where m'^ 
is zero, it is sufficient to obtain the required result that these contractions 
should have the property of uniformity, i.e. that the should be constants, 
over the region in question. 
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The proof the circularity theorem to this effect uses the fact that the 
Kilhng equations (4.1) and (4.2) imply corresponding higher order conditions 

V,V"P = -R^'X , VuVm^" = -R^'^m" , (4.20) 

from which, with the aid of (4.3) one gets 

V^icu^m,) = \e^,p„k''mPR\k^ , = le^,p^m''kPR\k^ . (4.21) 

In the absence of electromagnetic source contributions the material circular- 
ity conditions (4.17) alone are sufficient to ensure that the right hand sides 
of the foregoing pair of equations vanishes, which is evidently sufficient to 
establish the required uniformity. To show that the result remains valid in 
the presence of electromagnetic effects requires a little more work, starting 
from the group invariance conditions 

k'V.Fp^ = 2F,[pV^]k' , m^V.F,, = 2F,[,V,]m'^ . (4.22) 

These conditions, together with the Maxwellian field equations ( ?. ?) imply 
that the field vectors 

= F^,k^ , B^' = '-e^'P'^k^Fp, , (4.23) 

will satisfy 

V,{E,m'') = , V^(5^m,) = Ane^^p^k'^m^r , (4.24) 

in which not only the first but also the second of the right hand sides will 
obviously vanish wherever the current singularity condition (4.16) is satisfied, 
with the implication that E^m'^ and B^rrip^ will also both be uniform and 
therefore vanish 

E^nf = , B^'rUf, = (4.25) 

since they both obviously must vanish on the axis. The conditions (4.25) can 
appropriately be described as field circularity conditions, since they are suffi- 
cient for the corresponding electromagnetic contribution to the gravitational 
source to satisfy the analogue of the material circularity condition (4.17) so 
that its effect does not invalidate the conclusion that uj^mi, and ^l]pk'' will 
vanish also. 

The orthogonal transitivity property that is established in this way 
means that it will be possible, except where t<7 = 0, to choose the spacetime 
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coordinates in uch a way as to expess the metric in the standard Papeptrou 
form : 

ds^ = gijdx'dx^ + X{{d(f) - uodtf - w'^dt^], (4.26) 

where the coefficients X, uo w (as defined by (4.5) and (4.6) are indepen- 
dent of the "ignorable" coordinates ip and t, but functions of the two other 
coordinates x*, i = 1,2 whose locci of constancy are the orthogonal two 
surfaces predicted by the theorem. It can be seen that under these condi- 
tions the zamo trajectories themselves are orthogonal to the hypersurfaces 
on which t is constant, this hypersurface orthogonality condition (which, like 
(4.26) itself, would fail if the circularity conditions (4.16) and (4.17) were 
not satisfied) is what is meant by the statment that the zamos congruence is 
irrotational. 

To understand what happens on the zamosurface where the zamo 
worldlines become null so that w vanishes and the metric form (4.22) be- 
comes singular, we use the fact that the Frobenius conditions (4.19) on which 
it depends imply 

2w^[uVp]W^r = '^gt^Xp'^ut] ■ (4.27) 

This gives an equation for the gradient of orthogonal to the Killing vector 
surfaces of transitivity that is somewhat analogous to the one obtained for 
the gradient of the zamo angular velocity u directly from its definition, i.e. 

zi7[^^V^]CC7^ = rv'^Vi^War] , X^^iar'^pi^^ = 2ti7^m[^ V^m^] , (4.28) 

the noteworthy thing about both these equations being that their right hand 
sides vanish on the zamosurface where w"^ = 0. Since both w and uj are 
invariant under the group action their gradients everywhere must be orthog- 
onal to the surfaces of transitivity generated by the Killing vectors, whereas 
according to (4.28) they must actually be aligned with these surfaces of tran- 
sitivity on the zamosphere, conditions which can only be reconciled if they 
are both aligned with the unique combination of Killing vectors that is null 
on the zamosphere, i.e. with the zamo direction itself, and hence that they 
are aligned with each other. Provided the gradient V^CiJ^ is non zero, and 
so defines the direction normal it the zamosphere it obviously follows that 
(as can be shown with a little more care to be true in any case''^^^) the the 
zamosurface is a null hypersurface and that the zamo angular velocity uu 
has a uniform value say on it. This is the result that I refer to as the 
weak rigidity theorem}-'^^^ (the qualification weak being because it is bases on 
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a line of argument assuming axisymetry at the outset, in contrast with the 
stronger rigidity theorem of Hawkingt'^'^1'''] based only on an assumption of 
analyticity) . 

The uniform angular momentum value whose existence is thus estab- 
lished can be extrapolted off the zamosurface to a uniform value throughout 
space, 

V^n« = , (4.29) 
in terms of which can construct a unique Killing vector combination 

F = k'' + n^m^" , V(^4) = (4.30) 

which is characterised by the property of becoming aligned with the zamo 
direction where this direction becomes null, i.e. on the zamosurface which we 
now know to be itself null. This shows that (unlike an ergosurface which is 
typically timelike, and subject of course to the postuate that the circularity 
conditions (4.16) and (4.17) are satisfied) the zamosurface is automatically 
what I have called a Killing horizon^^^\ i.e. a null hypersuface whose null 
generator conicides with the generator of an isometry. 

Before going on to consider the global question of the identification of 
the locally defined zamosurface Killing horizon with the globally defined black 
hole event horizon, there are some further local properties of Killing horizons 
that can logically be derived at this stage. To start with it is apparent from 
the Penrose Raychaudhuri equation (3.11) for the null generator that since 
all the other terms vanish we must also have 

R^^Ft = (4.31) 

which subject to the material energy positivity postulate means that we must 
separately have 

= , T^^£/. = (4.32) 
where the electromagnetic part is given by 

(4.33) 

where a dagger symbol is used to distinguish quantities defined with respect 
to the corotating Killing vector field (4.30) from their analogues as defined 
with respect to the ordinary (asymptotically timelike) stationarity Killing 
vector (4.1) ( a distinction that is not necessary in the static case for which 
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they both coincide). If the material contribution is simply of perfect fluid type 
(2.25) subject to the inequalities p > , P > the first of the conditions 
(4.28) can be seen to give 

p + p = =^ p = o,P = 0, (4.34) 

i.e. there must be a vacuum at the horizon, while since both and are 
both by construction orthogonal to they cannot be timelike on the Killing 
horizon, so the reconciliation of (4.32) with (4.33) requires that they both be 
null there and hence proportional to the generator itself, i.e. 

E\/,] = , B\/,] = , (4.35) 

of which the first tells ust^^^ that the horizon is like a conductor whose equi- 
librium requires uniformity of the corresponding potential, an analogy that, 
since it was first noticed, has been developped in considerable detail'^^t^°°"^°^^. 
Explicitly we have 

4 = V^$^ , = A^F (4.36) 

everywhere, with the potential necessarily uniform over the horizon. For- 
mally, for any tangent vector to the Killing horizon we have 

= e'V^$t = . (4.37) 

This uniformity property of the angular velocity and of the po- 
tential <l>^ on a Killing horizons are prototypes'^^^t^^^t^^l for a less intuitively 
obvious uniformity'^'^"^''"^^"^^' , namely that of the decay parameter k whose 
definition by the general formula (3.7) is unambiguous now that the normal- 
isation of i'^ is fixed by (4.30). The fact that i'^ must satisfy the Frobenius 
orthogonality condition (3.2) on the horizon can be seen to mean that there 
must exist some vector on the horizon such that 

VJ^ = 2g[,£^i , Tg, = (4.38) 

where is an ingoing nullvector as introduced in (3.5). It is easy to see from 
the expression (4.38) that for any vectors t]^ lying in the horizon we shall 
have 

C^e, = , r]^i, = er^'^VJ, = ^ i^r^'^V.i, = , (4.39) 
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which shows that the Kilhng horizon is extrinsically flat (geodesically gen- 
erated) since it shows that rj^VuC,^ will automatically be tangenial to the 
horizon. A further derivation now leads (using i^^p^'^ = ^'^qu^p) to 

e^r/'^V.V/. = -{V ,inuWV pT,'' + r^'^V^rhoi^) = (4.40) 

Since the Killing vector property (4.30) by itself implies 

VpV^4 = i?^.p"4 , (4.41) 

we end up with 

^M-prnV = 0, ^ W^R,uprte = Q. (4.42) 

where 'g^^ is the projection tensor given by ((3.6), whose substitution then 
gives 

R^pte = R^.pJ'^ri'C ■ (4.43) 
Since the defintion (3.7) is clearly equivalent to 

K = -FrVJ^ = tq, [AAA) 

direct differentiation gives 

= -e (rrv,v,£^ + r(Vpr) v,£^ + «:£^v,r) (4.45) 

in which the last two terms cancel since they are respectively equal and 
opposite to i^C^Qu, so using (4.43) one finally obtains the simple result 

rV.K = -R,,rC ■ (4.46) 

The reasonning up to this point has been purely kinematic. If we now in- 
voke the Einstein equations we immediately obtain the required uniformity 
condition that for an arbitrary tangent vector xi^ to the horizon 

^"^V^K = (4.47) 

in the pure vacuum case. This "zeroth law of black hole mechanics" result t^^' 
can easily be seen to remain valid for a source free Einstin Maxwell vacuum'^^^t^^J 
since in that case, although R^^i, will not be zero Rfj^yd^ will be be proportional 
to the null tangent covector £^ which is all that is required to get (4.47) from 
(4.46). 



49 



5. Rotating Equilibrium States : the global problem. 

The kind of stationary equilibrium state towards which, subject to 
the cosmic censorship hupothesis, it seems reasonable to suppose an isolated 
gravitationally collapsing system would evolve, and that we shall understand 
to be meant by the qualification "well behaved black hole equilibrium state" 
is a stationary spacetime whose domain of outer communications (D.O.C.) 
is bounded to the future by a well behaved black hole event horizon. 

In the concrete example of the Schwarzschild and Riessner Nordstrom 
solutions we have seen that the D.O.C. is also bounded to the past by a 
well behaved "white hole horizon", but the latter corresponds to nothing 
that exists in a dynamical system that collapses from well behaved initial 
condition (being merely an artefact of analytic extrapolation to the past) 
whereas the black hole horizon in the stationary state really does correspond 
to the limit of the black hole horizon of the dynamically collapsing state. This 
is why only the latter can appropriately be postulated to exist as a defining 
characteristic of a well defined black hole equilibrium state, the existence of 
any other horizon in the past being something to be proven (if it can be) 
subsequently, but not to be postulated in advance. If it is already known 
in some particular case not only that both past and future event horizons 
exist but also that they have a well behaved Kruskal type crossover on a 
spacelike two surface then one can prove a result such as the "zeroth law" 
obtained at the end of Section 4 by a much shorter argument than was given 
there (see the accompanying lectures of Wald) but it is important for our 
present purpose to have shown that the result (i.e. the uniformity of k) can 
be established independently of any such assumption since it is needed as an 
an intermediate step in the line of reasonning that ultimately shows (at least 
in the generic vaccuum case, the more general question remaining open) that 
the assumed crossover really does occur. 

In order to have the right to utilise the local properties established 
for a zamosurface Killing horizon in the previous section we must show that 
the globally defined black hole event horizon forming the future boundary 
of the D.O.C. really is of this type. As a step towards this it, is convenient 
to use a 

lemma[i2l[ioi[57] 

giving a pseudo local characterisation of the D.O.C. 
whose original global definition is expressible for a stationary state as the 
specification that it consists of the intersection of the past and the future of 
the outer region where the trajectories of the stationary generator, d/dt = 
k^d/dx^ are actually timelike. The pseudo local characterisation is based on 
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consideration of where these generators are globally bradyonic, using this 
term to qualify any curve in spacetime with the property that any point x 
within it determines corresponding points x+ and x_ such that the part of 
the curve preceeding X- lies entirely in the past of x and the part of the curve 
subsequent to x+ lies entirely in the future of x. This condition is satisfied 
trivially by an ordinary timelike curve (for which a;_ and a;+ can be identified 
with the original point x itself) but it also includes a curve which, though 
locally spacelike, turns back towards itself in such a way that its evolution is 
effectively timelike in the long run. 

It is not difficult to see'^^"^''"^''' that although there may be an ergore- 
gion in which the Killing vector ceases to be timelike, its trajectories must 
always remain globally bradyonic throughout the D.O.C. It is even more ob- 
vious that in any connected region such that the trajectories there are all 
globally bradyonic, any one of them can be conected to any other by a time- 
like line with either orientation. This leads to a /emma^^^'I^^lf^^l characterising 
the D.O.C. as the maximal connected extension of the outer region where the 
stationary trajectories generated by /c^ are timelike such that the stationary 
generators remain at least globally bradyonic. As an immediate corollary it 
follows that can never be timelike on the boundary (including the black 
hole horizon) of the D.O.C, while on the other hand the stationarity gener- 
ator k^^ can never vanish nor have any closed trajectory within the D.O.C. 
The latter conclusin means that k^ must always be linearly independent of 
m'^ in the D.O.C. and hence that within the D.O.C. the Killing bivector w'"^ 
can never be degenerate except on the axis = 0, which means that ro^ can 
vanish only on the axis or where the bivector, and thus also the corresponding 
zamo direction, becomes null. 

Let us now designate by Z the maximal connected extension of the 
outer region where the stationary trajectories generated by k'^ are timelike 
such that the local condition that the zamo trajectories be timelike within it 
is satisfied. Our aim is to show, subject to very weak assumptions that this 
locally defined domain will be identifiable with the globally defined D.O.C. 
What is obvious is that, since any connecting curve within Z will be contin- 
uously deformable into a timelike curve (by the group action along the zamo 
direction at each point) Z must certainly lie entirely within the D.O.C. The 
definition of Z means, according to (4.10), that Z is characterised by w"^ > 
except on the rotation axis where w'^ = and that we must have zu"^ = 
everywhere on the boundary Z of Z. By the conclusion of the previous para- 
graph, this implies that except on the axis Z lies on the locus where the 
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zamo direction is null, and hence by the results of the previous section that 
provided the circularity postulate is satisfied ; the connected components of 
Z must be Killing horizon and thus null hypersurface segments, each, by 
continuity, with uniform time orientation. 

To complete the demonstration that Z can be identified with the 
D.O.C, we must now invoke the further postulate that the latter be sim- 
ply connected, which means that a (hypothetical) connected component V 
say of the complement of Z in the D.O.C. must have a boundary T> within 
the D.O.C. that is itself connected. It follows that T> would have to be a 
null hypersurface segment with uniform time orientation which means that 
V could be reached from Z only by future directed timelike lines, or only by 
past directed ones, but not by both kinds as would be required for V to lie 
within the D.O.C. so in order to avoid a contradiction it must be concluded 
that V is empty, and thus that Z covers the whole of the D.O.C. as required. 

The foregoing demonstration to the effect that the Killing bivector 
wf"' is timelike throughout the D.O.C. means that the two surfaces orthogo- 
nal to the Killing vectors that were shown to exist by the circularity theorem 
of the previous section will correspondingly be strictly spacelike there. More- 
over the reinforcemnt of the circularity postulate introduced in Section 4 
section by the simple connectivity postulate introduced in the previous para- 
graph implies that these orthogonal two surfaces will be constructible not 
just locally but globally. Since any such surface differs from flatness only by 
a locally variable conformal factor, S say, it follows that the space coordi- 
nates in the general Papaptrou form (4.26) may be chosen more specifically 
to be cylindrical type coordinates p, z say, in such a way that the metric will 
be expressible in the form 

ds"^ = T.{dp^ + dz^) + X{{d<i) - udtf - w^dt^} (5.1) 

which will be globally valid over the entire D.O.C. except for the familiar 
degeneracy on the axis where X and w"^ vanish, their vakues elsewhere being 
strictly positive, as is the value of E everywhere. An analogous form is ob- 
tainable for the vector potential using the invariance conditions (4.22) which 
imply the existence locally, and hence by the simple connectivity postulate 
globally, of scalars $ and B such that 

Fp.k'' = Vp$ , Fp^m'^ = Vp5 , (5.2) 

in terms of which, using the consequence (4.25) of the circularity postulate. 
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it can be verified that the gauge may be chosen in such a way that 

ApdxP = Mt + Bd^ . (5.3) 

It is of course to be understood here that the new coefficients S $, 5, hke 
X, u, and w as introduced previously, areall functions only of p and z only, 
i.e. the stationarity and axisymety is mad manifest by their independence of 
the ignorable coordinates (f) and t. 

Up to his stage the analysis has been sufficiently general to cover a 
wide range of conceivable black hole configurations with external matter rings 
for which explicit analytic solutions are not available, but from this point on 
we shall restrict our attention to the globally source free case for which it has 
long been knownt'^^l['^^l[^°l['^^l['^^l that the Kerr Newman class of solutions pro- 
vide explicit examples. The purpose of the systematic step by step approach 
whose development I have been describing is to solve the problem of whether 
there can exist any others. It will be shown below that subject to the pre- 
ceeding assumptions, including notably that of a simply connected topology 
for the D.O.C, it can be shown, by an argument that has been able to be 
made completely watertight only comparitively recentlyl^^^I^^lt^^l , that these 
known solutions are indeed the only ones for the strictly source free Ein- 
stein Maxwell equations. However the problem remains wide opent^^l for the 
slightly more general case of solutions of the source free Einstein Maxwell 
equations with cosmological A term : such solutions cannot of course be 
asymptotically flat, but I have discovered[''^]t''^]t^°] a wide class of asymptoti- 
cally De Sitter solutions (one of the first cases for which the CP. technique 
described in Section 2 proved quite indispensible for providing an under- 
standable global description t^'^l'^^l. The problem that remains unsolved is the 
extent to which these known asymptotically De Sitter black rotating black 
be whole solutions are unique. 

The reason why the results that follow have not yet been generalised 
to allow for a cosmological A term in the generalised source free Einstein 
Maxwell system 

R^. - ^Rg^'u = 87rT/„M + Ag^, (5.4) 

is that the next step uses the trace, not over the full four dimensional system 
(5.4) but over its restriction to the two dimensional surface of transitivity 
generated by the Kiling vectors for which due to the cicularity condition 
te electromagnetic contribution cancels out so that substitution of the form 
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(5.1) simply gives 

—V^w = 2A , (5.5) 
azu 

where is the Laplacian that is defined with respect to of the flat two 
dimensional metric dp'^ + dz'^. 

In the absence of the cosmological term, (5.5) tells us the is a 
harmonic function on the conformally fiat spacelike two surfaces which means 
that using the freedom to make conformal adjustments one can choose the 
coordinate system so as to identify it with p : 

A = ^ w = p . (5.6) 

In such a coordinate system (5.1) reduces to the specialised Papapetrou form 

ds^ = J:{dm^ + dz^) + X{{d<p - ujdtf - w^df} . (5.7) 

On the other hand the presence of a A term suffices to block the apparently 
innocent but actually crucial step from (5.1) to (5.7), without which none of 
the work that follows can be carried through, the discovery of an alternative 
route being thus left as a challenge for future work. 

As has been well known since early work^'^^lt^'^l on stationary axisym- 
metric systems in other contexts, after the variable to has been thus taken out 
of the list of unknown variables by its promotion to the status of a "known" 
coordinate variable, the system of source free field equations reduces to a 
decoupled system just for the two metric variables X and u together with 
the two electromagnetic potentials $ and B, together with a separate equa- 
tion that can be solved afterwards to obtain the remaining variable, i.e. the 
conformal factor S by a direct quadrature, with the constant of integraton 
fixed by the condition S — * 1 at large asymptotic distance. This means that 
the four dimensional black hole equilibrium problem from which we started 
reduces now to a two dimensional boundary problem for the fields X, u, $, 
B as functions of the independent variables w and z. 

The conclusion that (5.7) is not just valid locally but covers the 
entire D.O.C. with w (ranging from to oo) and z (ranging from — oo to 
oo) as globally well behaved cylindrical type coordinates, depends on the 
simple connectivity postulate and the use of a specialisationt^^l of Morse 
theory to exclude the possibility of critical points of w. Ordinary Morse 
theory establishes that the number of maxima plus the number of minima 
minus the number of minimaxes is a topological invariant subject to fixed 
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boundary conditions and the assumption that no degenerate critical points 
occur. In the harmonic case it is possible to make a stronger statement since 
maxima and minima cannot occur while degenerate critical points do not 
need to be assumed to be absent since their presence can easily be allowed 
for by labelling them with an apropriate positive degeneracy index. The 
resulting theorem^^^l states that the index weughted sum over all critical 
points including possible degenerate ones is a topological invariant which in 
the present application can be seen (by considering any special case such 
as the Scchwarzschild solution) to be zero. Since the index is always of the 
same sign in the harmonic case no cancellation is possible, so the fact that 
the total is zero in the case under consideration makes it possible to deduce 
with certainty that there are no critical points at all, degenerate or otherwise. 

Having got to this point we introduce a further topological simplifica- 
tion postulate to the effect that we are dealing with only a single topologically 
spherical black hole. Our simple connectivity postulate has already excluded 
conceivable toroidal black hole configurations but has left open the possibility 
of having several topologically spherical black holes lined up on a common 
rotation axis. A certain amount of work has been carried out, particularly by 
Hawkingt^'^] and Gibbons'^^' towards showing that such configurations are im- 
possible except in the extreme Papapetrou Majumdar limit t^^lt^'^^t^^^l^^l'^^ll®^] of 
maximally charged non rotating configurations in which electromagnetic re- 
pulsion balances gravitational attraction, but we shall not go into the study 
of such exotic topological possibilities here. Assuming then that the black 
hole topology is of simple spherical type we can fix the cylindrical coordinate 
system symmetrically with respect to the black hole by taking the poles at 
which the horizon meets the rotation axis to be given by opposite values of z 
which, it can easily be seen, must be given explicitly by in terms of the area 
and decay constant of the horizon by 



Leaving aside the awkward and still only partially understood spe- 
cial case for which the horizon is degenerate in the sense of having k = 
(corresponding to zero temperature in the thermodynamic limit) we can con- 
veniently proceed by replacing the cylindrical coordinates t<7, z by ellipsoidal 
type coordinates A, according to the specifications 



^ = ±c , 






(5.9) 
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which are such as to arrange that the horizon is now given by the hmit X c 
while the two disconnected ("north" and "south") parts of the symmetry 
axis in the D.O.C. aregiven respectively by /i — > ±1, the whole D.O.C. being 
covered by the coordinate range c<A<oo, — l</i</i. 
In this system the metric takes the form 

ds^ = Eds^ + + 2Wd^dt - Vdf (5.10) 

for a conformally flat space metric given by 

^ 2 dX^ dfi"^ 

with 



=Tf^ + T^ (5-11) 



H = (A2-cV)S W = Xu;, V = X-\w^ -W^) . (5.12) 

In terms of the two dimensional covariant differentiation operator V defined 
in terms of the known conformally flat metric (5.11), the system of indepen- 
dent source free Einstein Maxwell equations reduces to a pair of Maxwell 
equations 

v| — (V$ -cjVfi)) = , (5.13) 



v|^VB + -(V$-cjV5)1 = , (5.14) 
together with a pair of Einstein equations 

— Vuj + — {V<^ -uWB)\ =0 (5.15) 
zu zu J 

V — VX + J ^ + — |V$ - iuVB\^ + —\VB\^ = . (5.16) 

La } to to X 

Although this system is singular on the axis where X and to both vanish, 

and also on the horizon where to also vanishes, it is garanteed to be regular 

everywhere within the half plane under consideration where we have 

to>0 , X >0 (5.17) 

the latter inequality being derived from the causality postulate. This is the 
motivaion for having used a formulation giving the leading role to tha ax- 
isymmetry Killing vector m'^ rather than using the more traditional approach 
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giving the leading role to the stationarity Killing vector k'^, which would have 
given an analogous system but with V turning up instead of X in the denom- 
inators, which would have the seriously inconvenient consequence of making 
the system singular on the ergosurface that generically occurs within the 
D.O.C. 

The foregoing system can be made more tractable by performing the 
analogue of the transformation introduced originally for the traditional for- 
mulation based on k'^ rather than by Ernst t^"'. This is done by first using 
the Maxwellian equation (5.13) to justify the introduction of a stream func- 
tion type electric potential E given by 

^ U - " aA J = - ^ ) a;; ' ^ - j = - ' 

(5.18) 

and by using the Einstein equation (5.15) to justify the introduction of an 
analogous rotation potential Y given by 

..2duJ , nAdY SB dE] 

^duj f dY dB dE] 

Using the new potentals E and Y to replace $ and u one obtains the Maxwell 
equations in the form 

V {^Vsj + ^{(VF + 2EVB - 2BVE} ■VE = , 

V {^Ve} - ^^('^^ + ^^^^ - ^BVE} . VB = , (5.20) 
V {^(VV' + ^EVB - 2BVE)^ = , 
V {^Vxj + ^ {|VX|' + |Vr + 2EVB - BVE\^} + 

^{|V^|' + |V5p} = . (5.21) 

In terms of this new system the asymptotic boundary conditions for 
regularity at large distance, i.e. as A — > oo, are more complicated than in the 
traditional approach, being obtainable'^"] as 

\-^X = (1 - + o(A-^) , E = -Qfi + 0(A-^), 
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r = 2J/i(3 - /i') + 0(A"^) , B = 0{X^^). (5.22) 

where J is the asymptoticaly measured angular momentum about the rota- 
tion axis while Q is the total charge, and where the requirement that the 
total magnetic monopole should vanish has been taken into account. The 
asymptotic mass M does not appear explicitly, but it is implicitly fixed by 
the overall scale which is determined by the choice of the parameter c. In 
compensation for this rather inhabitual degree of complication in the famil- 
iar large distance limit, we get extremely simple boundary conditions in the 
less familiar limit at the horizon as A — > c, the only condition here being 
that the unknowns E, B, X, Y should be regular as differentiable functions 
of the ellipsoidal type coordinates A and /i. The most mathematically deli- 
cate boundary conditions (for which however no physical considerations or 
parameter values are involved) are those for geometrical regularity on the 
rotation axis fi = ±1, i.e. for the limit (1 — mu'^) 0, which are given by 

X = 0(1-/.^), ^-^^-^^ = ^ + 0{l-^,^). (5.23) 

After I first obtained this system I succeeded deriving a pure vacuum 
"no hair theorem" '^^^It^'^l which Robinson was able soon after to generalise 
to the full electromagnetic case^^^^, establishing that the solutions belong 
to descrete families each depending continuously only on the three relevant 
physical parameters involved in the boundary conditions namely Q, J, and 
C, of which the latter fixes the overall scale and thus implicitly the mass 
M. Our method was to equate a certain divergence to a positive definite 
function of the infinitesimal difference between nearby solutions for the same 
parameter values and hence (using the boundary conditions) to show that 
the latter must vanish. One such family consisted of course of the already 
known Kerr solutions (subject to the condition > + J'^/M'^) and in 
view of various restrictions on special limits such as that of spherical symetry 
it seemed unlikely from the outset that any others existed. Nevertheless it 
was necessary to wait several years before it was established beyond doubt 
that they do not. 
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The way that Robinson and I had constructed the divergence with 
the miraculously positive definite form we needed was based on a purely trial 
and error approach whose success in the electromagnetic caset^''' required 
a veritable algebraic tour de force. Robinson even succeeded in using the 
trial and error metheod to construct a finite difference generalisation'^^] that 
complelty established the uniqueness of the original Kerr black hole solutions, 
with Q = 0, Af2 > J2/m2, for the pure vacuuum case, but the trial and 
error method was never able to cope with the finite difference case in its 
full electromanetic generality, and so the complete solution had to await the 
introduction of new and more sophisticated techniques by Bunting'^^lps] g^j^^j 
Mazurt^^l. The Bunting method is of great interest in its own right, being 
potentially useful for much more general problems'^^' . I shall however restrict 
myself here to the description of the Mazur method which is more specialised 
but more explicit. 

As soon as I obtained the system given above I noticed that the field 
equations (5.21) and (5.22) have the striking feature'^"' (which does not apply 
to the traditional system defined in terms of k'^ rather than m^) of being 
derivable from a Lagrangian integral 

X= J CdXdfx (5.24) 

that is positive definite with the comparitively simple form 

^ |VX|2 + |VF + 2EV5-25VE|2 |VE|2 + |VfiP 

C = ^^ '- + 2 ' (5.25) 

but neither I nor Robinson had seen how to exploit this directly. The break- 
through by Mazur was based on work by Geroch'^^^ and Kinnersley'^'^^ who 
showed that Ernst type systems can be interpreted as belonging to a class 
of non-linear a models whose field equations are equivalent to a partially 
redundant set of ordinary divergence type conservation laws of the form 

VJ = (5.26) 

where J is a matrix vector constructed according to the prescription 

J = tu~^<P^^ ■ W (5.27) 

where w is the known positive weight function given by (5.9) (which does 
not appear in traditional a models but whose appeareance here adds no sig- 
nificant complication) and where ^ is a hermitian matrix function with the 
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important property of being positive definite in the present case, its compo- 
nents being given by 

$dfe = Vdb + 2vaVb (5.28) 
(using a bar to denote complex conjugation and placing a dot on conjugately 
transforming indices) where T]ab is just the the fixed Minkowski type hermitian 
metric in diagonam form with signature (—1,1,1) for the space of three 
dimensional complex vectors, in which the field vector Va is given in terms of 
the complex Ernst type variables 

e = -X + iY - ij^j , ij = E + iB (5.29) 

by 

{v„ v„v,) = -1,6 + 1,2^) (5.30) 

which is such as to make Va automatically "timelike" in hermitian space, with 
unit normalisation given by 

V^'vaVf, = -1 , (5.31) 

which is sufficient to guarantee the required positivity of the hermitian matrix 
given by (5.28). 

The preservation of the field equations by the SU(2,1) action (leaving 
rj^jj invariant) that Kinnersley thus made manifest can also be seen to extend 
"off shell" in the sense that our Lagrangian (15.25) is also invariant, as can 
be seen by rewriting it as 

C = 2\r]^%aVv,\^ - r]''\Vva)Vv, = i^.,tr{J* ■ J^} (5.32) 

where Qij {i,j = 1, 2) is the positive definite two dimensional metric given by 
(5.11) that is used for the specification of V. 

The Mazur method of establishing the uniqueness of the solutions of 
such a system, subject to appropriate boundary conditions such as are given 
in the present case, is essentally dependent on the positive definiteness of 
both <P and Qij. The objective is to prove the vanishing of the difference 

i = - <P[o] (5.33) 

between any pair of matrices representing conceivably distinct solutions. Te 
vanishing of this difference is evidentlyequivalent to the vanishing of what! 
refer to as the deviation matrix, 

A = ^>[i] ■ - 1 = i ■ (5.34) 
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where 1 is the unit matrix (in the complex 3 space). The gradient of this 
deviation matrix will evidently be given by 

wVA = <P[^ ■ J ■ <l>py 1 , J = J[i] - J[o] . (5.35) 
Taking the difference we obtain 

V(tuVZ\) = ^>[i] • |V J + t^-^ft,(Jii] ■ Jfi] - 2J[i] ■ Jfo] + Jio] ■ J[o])} ■ ^[oT' 

(5.36) 

The next step is to use the hermiticity property 

<p = <p* =^ J* = <P-J-<p-^ (5.37) 

(where the asterisk denotes the complex conjugate of the transpose) to rewrite 
the quadratic terms in (5.36) as 

On taking the trace of (5.36) we thus obtain the scalar identity 

V {wVtT{A}) - tr{<?[oyi ■ ^>[i] ■ J} = w~'g,, ti{^^^' -3*^- ^>[i] ■ • 

(5.39) 

This Mazur identity includes as special cases the identities found for the 
linearised or uncharged limits by Robinson and myself '^^^ t^^^ t^^^ t^^^ using a less 
systematic approach. 

It follows directly from the form (15.26) of the field equations that the 
current difference satisfies ^ 

V J = (5.40) 

and hence that the left hand side of the identity (5.39) will reduce to a 
divergennce whose integral can be converted using Green's theorem to a 
surface contribution which will vanish, 

^ dSi wg'^VjitiA) ^ (5.41) 

subject to the appropriate boundary conditions which can be verified (using 
particular care for the axis where diverges) to be in fact satisfied in 
this case. Under these conditions one can deduce that the right hand side of 
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(5.39) vanishes since (by the positivity of ^) it is clearly a positive definite 

± 

function of J which must therefore vanish, i.e. using (5.35) we get 

J = ^ A = C (5.42) 

where C is some constant matrix. Since the boundary conditions as A ^ oo 
ensure that Z\ — > there one ends by getting 

C = ^ ^>[i] = ^>[o] (5.43) 

which finally establishes the required uniqueness. 

6. Special Properties of the Kerr Newman Vacuum Solutions. 

The theorem obtained at the end of Section 5 establishes conclusively 
that in the source free electrovac case there are no (topologically simple) 
stationary axisymmetric asymptotically fiat black hole solutions with non 
degenerate {k > 0) horizon other than those of the Kerr Newman family'^^^l 
as restricted by the condition 

> , = M'^ - a'^ - Q\ a = J/M , (6.1) 

where this parameter c represents the value for these solutions of the quantity 
introduced more generally by (5.8) . This includes, for Q = 0, the pure 
vaccuum family of Kerr'^'^^l solutions, whose black hole nature, in the allowed 
parameter range > a^, was first clearly recognised by Boyer[''^][''^]. (For 
the degenerate limit for which c = or equivalently k = 0, corresponding to 
a horizon at zero temperature in the thermodynamic analogy, the problem 
has still not been completely solved, but it is known that in the non rotating 
case J = there is a classt^^^^^^'f^^J^^^lf^^J^^^l of electrically balanced solutions 
with = that is much more general than the corresponding Reisner 
Nordstrom subset within the Kerr Newman familly.) Taking full advantage of 
the very special properties that will be briefiy surveyed below, a long series 
of investigations, of which the first was that of Vishveshwarat^^^' and the 
most recent that of Whitingt^^^ confirm that the Kerr solutions specified by 
(6.1) are effectively stable against the all the most obviously relevant kinds 
of perturbation. 

In terms of ordinary coordinates r, 6 introduced by 

\ = r- M , fi = cose (6.2) 
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the explicit solutions for tha Ernst type variables are given by 



X = +a + \ sm 6 

f 2.x sin2^[gV - Ma^sin^^l 1 
r = <^ M 2 + sin'^ ^1 i 2a cos ^ 

Q(r2 + a2)cos 6 -Qra sin^O 

-C/ = 7: 7: — , -D = — — - . (D-'jj 

+ a^cos^t' r"^ + a^cos^y 

Going back to ordinary metric and electromagnetic potential components 
gives the rather simpler forms 

y^^_ 2Mr - ^ Qr 



+ a^cos^6' ' + a^cos^6' ' 

^ = — — 5 -T^^ ' ^ = + a^cos^e . (6.4 

+ a^cos^y 

Subject to (6.1) these solutions do in fact have turn out[^''l['^'^l['^^l[^°l to 
have the property (which was not assumed in advance in the approach out- 
lined above) of having a well behaved Kruskal type horizon crossover when 
analytically extended towards the past. When analytically extended to the 
interior they exhibit many amusing but one presumes physically irrelevant 
features such as a time machine [''^^t^^l in the region beyond the Cauchy hori- 
zon, which (in this more general case, as in the Reissner Nordstrom case 
discussed in Section 2 and in more detail in the accompanuing lectures of 
Israel) is a sign of instability occurring at r = r_ , using the standard abrevi- 
ation 

r± = M±c. (6.5) 

In these solutions the three quantities whose uniformity over the black 
hole event horizon at r = r+ was guaranteed in advance by the results of 
Section 4, namely the decay constant k, the limiting value of the zamo 
angular velocity u, and the value $^ of the comoving potential $''" = $ + 
Q"B, will be given respectively by 

2Mr+-g2' 2Mr+-Q2' ^ 2Mr+ - • ^'^""^ 
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Neither the physically motivated global geometrical approach that we 
have followed here, nor the even the analytical approach (based on the as- 
sumption of a special form for the Weyl tensor) used originally by Kerr and 
Newman makes it at all obvious in advance that the solutions should have 
such remarkable special algebraic properties as they actually do. The special 
simplicity of the metric that is finally obtained can be made most directly 
manifest by expressing it in tetrad form 



terms of a preferred canonical separable tetrad which I found very early in the 
history of the study of the Kerr Newman solutions when the separablity of the 
Hamilton jacobi equation for free particle trajectories and of the scalar Klein 
Gordon wave equation was first brought to light '^^"^''1. However although it 
is useful generally[^°^l[^°^][^°^l[^°^l and is even more strongly to be preferred 
for obtaining separability of higher spin wave equationst^°^"^^°lt^^^"^^^' , this 
canonical tetrad has unfortunately tended to have been neglected by workers 
on higher spin separabilityt^^^l'^^^lt^^^l in favor of other tetrads (particularly 
that of Kinnersleyf^^^l which has advantages in other contexts but not in this 
one) thereby making many unavoidably heavy calculations'^^^^ even longer 
and more complicated than necessary. To get maximum algebraical symmetry 
it is necessary to replace the geometrically defined time and angle coordinates 
t and as introduced in the previous section (and which in the specific 
context of the Kerr solutions are commonly referred to as the coordinates 
of Boyer Linquist'^^' in order to distinguish them from the rather different 
time and angle coordinates used by Kerr himself '^"^^ ) by very closely related 
ignorable coordinates t and say, and to introduce a recalibrated angle 
coordinate q say in place of /i, according to the specifications 



ds'^ = dxf'dx'' (-€"„€' 



(6.7) 



t = t — a(j) 



= a 



1 



q = a cos 6 . 



(6.8) 



The canonical maximally symmetric tetrad is then expressible as 
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where and Ag are quadratic functions repectively of r only and of q only, 
the symmetry between these variables being broken only by the presence of 
a linear mass term in the former but not the latter : 

A^ = - 2Mr + , A, = - . (6.10) 

The electromagnetic field potential is expressible in the even simpler form 
(proportional to the first of these tetrad forms as 

A^dx^ = [^:^^^^ {q^d4> + di) (6.11) 

At the expense of violating the asymptotic boundary conditions the 
algebraic symmetry between r and q could be made complete'^''' by includ- 
ing appropriate gravimagnetic and electromagnetic monopole terms. The full 
symmetry is however already manifest in the corresponding expression'^^^ for 
the crucially important Killing Yano 2-form /^i, given by 

f^.udx^'dx" = qdr A {di + q^d(j)) + rdq A {di - r^d(j)) , (6.12) 

whose existence, as a solution of the Killing Yano equations 

Uu = 0, Vp/^. = V[p/h, (6.13) 

underlies the remarkable hidden symmetries of the Kerr Newman family and 
is by itself sufficient to characterise them completely among asymptotically 
fiat electrovac solutions. 

Why the purely local characterisation given by (6.13) should give the 
same result as the global boundary conditions for a black hole equilibrium 
problem remains mysterious, but once it is known that (as was first revealed 
by the work of Penrose and his collaborators [us] [us] [120] [121] [122] ^gjj^g spino- 
rial approach) there is a non zero solution of (6.13), most of the other special 
properties of the Kerr solutions can be obtained by more or less straightfore- 
ward deduction without the intervention of any further independent miracles, 
starting with the stationarity, whose generator is given'^^^^ directly by 

p = l^M^P'^V./,. , (6.14) 

which, by (6.13) will automatically satisfy the Killing equation (4.1). In ad- 
dition to this "primary" Killing vector (6.11) also'^^^lt^^^'^^^l ensures the exis- 
tence of a "secondary", generically independent, one given by 

= a^^k'' , = f%r, , (6.15) 
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which again satisfies Kilhng's equation as a further automatic consequence 
of (6.13) and which turns out when evaluated to be given as a hnear combi- 
nation of the "primary" , stationarity generating Kilhng vector k^^ and of the 
axisymetry generating Kilhng vector (as distinguished by having closed 
trajectories) that is interpretable'^^^^t^^^^ as generating rigid rotations with 
angular velocity Q = = M/J about the axis : explicitly it satisfies 

V(^/i^) = , hf" = a^k'' + am^ . (6.16) 

The pair of independent Killing vectors thus obtained will of course give rise 
to a corresponding pair of quantities k^Uy and h^Uy (or k^Uy and m'^Uy) that 
are conserved allong solutions of the geodesic equations (4.8). The "hidden 
symmetry" corresponding to the existence of the "fourth" constant of motion 
that is needed to provide a complete set of first integrals of the equations 
of motion (the third being given trivially just bu m'^m^ itself) is given by 
the tensor a^'^ as defined in (6.15) which is an ordinary symmetric Stackel 
Killing tensor in the sense that, again as as automatic consequence of (6.13), 
it satisfies the conditions 

a[f^u] = , V^pQ/^y) = (6.17) 

which evidently suffice to ensure that the quadratic combination a^yU^u'^ 
will indeed be constant allong solutions of (4.8), thus providing the "fourth" 
constant that is needed to make these equations completely integrablel^^^I^^^. 
Moreover this tensor a^'^ is not just a Killing Stackel tensor in the weak sense 
of satisfying (6.17) but is automatically, by the integrability conditionst^^^t^^^l 
for (6.13), is automatically a Killing tensor in the strong sense^^^^' that it also 
satisfies 

a^^^Ry^p = (6.18) 

which is the suppmementary condition needed in conjunction with (6.17) 
to ensure that the corresponding self adjoint differential operator commutes 
with Dalembertian (scalar) wave operator, i.e. 

[VX'V.,V^V,] = , (6.19) 

Just as the scalar Dalembertian can be thought of as a sort of square of the or- 
dinary first order dirac operator 'j^'V^ acting on 4-spinors, so analogously'^^^^t'^^^f^^^ 
the operator V^a'^'^V^ can be thought of as a sort of square of a first order 
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generalised spinor angular momentum operator L that commutes with the 
Dirac operator 

[L,j^V,] = , L = z7/.(7-V/V. - k^) , (6.20) 

in the same way as does the ordinary Kosman'^^''' energy operator K for any 
Killing vector, 

[K, ^^V^] = , K = zPV^ + jli^'^Y^K, zV^] . (6.21) 

Just as the commutation law (6.19) is interpretable as resulting from 
the separability of the scalar Klein Gordon equation'^^1 so analogously the 
commutation law (6.20) is interpretable as resulting from the separability of 
the Dirac equation'^^^^ . In addition to these cases of separability for massive 
particle wave equations, the Kerre Newman solutions are also characterised 
by analogous separability properties for massless higher spin wave equations, 
including notably those for the separate electromagnetic and electromagnetic 
perurbations that are relevant to stability analysist^^^''^^^^, though so far not 
for the case when the electromagnetic and gravitational perturbations are 
coupled as will be the case for the generic perturbation in the charged Kerr 
Newman case. The analysis of such wave equations is carried out most con- 
veniently by using not the orthonormal version but the corresponding null 
version of the canonical tetrad, the latter being given in terms of the former 
by a transformation of the standard form 

= + eV) ' 4 = 7f - eV) , = ^(e°^ + ze^) , (6.22) 

where and £^ are null and 2;^ is complex. Let us consider together the case 
of an ordinary complex scalar field say, the case of an ordinary Maxwell 
field Ffj_iy, for which we take a priviledged pair of complex tetrad components 

= , = F^^z^r , (6.23) 

(where for reasons of notational convenience that will become obvious we do 
not use the traditional counting system which would label these components 
as $(, and $2 )' ^^cl finally the case of a gravitational perturbations with 
Weyl tensor C^j^^'cr for which we similarly take the priviledged pair of complex 
tetrad components 

^2 = -C^^.J^'z-'e^z'' , = -C^,pJ>'z''iPz" , (6.24) 



67 



(which in the traditional counting system would be labeled and Then 
the upshot of the many studies referred to above is that^^'^^' corresponding 
field equations are separable, in terms of our present notational system, by 
setting 

(r - = X,(r)F,(g)e-*(^*-*'^) (6.25) 

(where the helicity index s runs over the values 0, ±1, ±2) with the resulting 
separated equations having the form 

i d ^ d (Er^ + is(r - M) - ^ ] 

where K is the separation constant whose existence expresses the hidden 
symmetry and where it can be seen from the relation 

Et-^^ = Et-^(p (6.27) 

that E is interpretable as the ordinary energy associated with the "primary" 
Killing vector while $, which is analogously associated with the secondary 
Killing vector h^^ is related to the ordinary angular momentum constant $ 
associated with the axial Killing vector by the simple relation 

$ = a$ - a'^E (6.28) 

In the case s = the above form agrees directly with what is obtained in 
the limit of vanishing perticle charge and mass, e = m = from my original 
separated form of the Klein Gordon equation, but in the higher spin cases 
it differs from the forms originally obtained by Teukolsky due to his use of 
a non canonical (Kinnersley'^^^' type) tetrad which lead to an unnecessarily 
complicated form in which the symmetries manifest in the version above ( 
the helicity symmetry between s and —s, and the almost perfect algebraic 
symmetry between r and iq) are all spoiled by subjecting (6.9) to a symmetry 
violating tetrad transformation consisting of a combined boost and rotation 
of the form 
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\/r^ + q 

The dazzling prestige conferred on the Kinnersley tetrad by its successful 
use in the original discovery of the higher spin separability has unfortunately 
blinded many workers to the fact that the separation works even more effi- 
ciently in terms of the original canonical tetrad, with the result that many 
published calculations^^'^ are at least twice as long as necessary The canon- 
ical tetrad has however tended to come back into use in more recent work, 
whose achievements include separation of the equations of parallel transport 
of a tetradt^°^"^''^]t^°^l and the equations for the stationary equilibrium of cos- 
mic stringsI^^^Jt^^^l'^'^'^] of of certain simple kinds (not just the ordinary Goto 
Nambu kind but also the more general non dispersive model allowing for the 
averaged effect of noise or wiggles) whose mechanics will be explained in the 
following sections. 

7. Basic Brane Mechanics. 

The purpose of the last part of this course is to give a brief introductory 
overview (and some illustrative applications in the context of cosmic strings) 
of the general principles of brane dynamics using a recently developed fully 
covariant approach'^^"^^' that avoids the use of excess mathematical bagage 
(such as the use of distribution theory and specially adapted coordinates for 
separate subsystems) that may be useful for detailed calculations in specific 
applications, but that would obscure the simplicity and generality of laws 
such as the general equation governing the extrinsic motion of any brane, 
which is expressible in the formalism set up below (using underlining to 
distinguish quantities defined with respect to a p-brane under consideration 
from any higher dimensional analogue that may also be relevant) in terms of 
its stress momentum energy tensor T^*^, and its second fundamental tensor 
Ky.y'' in the form 

T'^^i^^/ = 1" , (7.1) 

_L 

where f is the total orthogonally projected force contribution (such as that 
of the wind on a sail, or of an external electromagnetic field on the current 
in a cosmic string) from the various external systems (if any) with which the 
brane may interact. 

Following an increasingly popular usage'^^^''^^^^, the term brane is used 
here to designate a physical model of the category that includes continuous 
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media and point particles as extreme cases, with ordinary membranes (from 
which the term is derived) and strings as the only other possibilities in a 
4-dimensional background. Generally, a (p — 1) brane is to be understood 
to be a dynamical system defined in terms of fields with support confined 
to a p dimensional world sheet surface 5 in a background (fiat or curved) 
spacetime manifold of dimension n > p. The extreme case, with n = pis that 
of a continuous medium for which the confinement condition is redundant. 
The use of this concept makes it possible to give a unified description of basic 
properties that are common to a very wide range of physically diverse phe- 
nomena. A simple and very important example is the universal rule that (as 
a consequence of (7.1) and independently of the nature of any external forces 
so long as their coupling does not involve gradients of internal field variables) 
the condition for a (contravariant) vector 77^ say to be an extrinsic bicharac- 
teristic vector, i.e. to be tangent to the direction of "group" propagation of 
localised wave packets of small extrinsic displacements of the localisation of 
the world sheet (which of course is meaningful only for n < p) with a corre- 
sponding characteristic covector x^l normal to the direction of the associated 
"brane wave" sheets, will be givent^^^l simply by 

= T^^x, , T^^x.Xu = . (7.2) 

The hyperbolicity condition to the effect that the characteristic equation 
(7.2) should define a real characteristic cone provides a restriction (trivial 
for a point particle and reducing just to a requirement of positivity of the 
ordinary tension T in the case of a string'^^^^) that must always be satisfied 
as a condition for local stability except of course in the case p = n of a 
continuous medium for which there is no geometric possibility of extrinsic 
perturbations, which is why an ordinary perfect fiuid with postive pressure 
P can be stable after all, despite the fact that it is elliptic (with no real roots) 
as far as the criterion (7.2) is concerned. 

In an ordinary spacetime with n = 4 a continuous medium (with 
p = 4) counts as a 3-brane, the other possibilities being that of a membrane 
model (with p = 3) which counts as a 2-brane, a string model (with p=2) 
which counts as a 1-brane, and finally at the other extreme, a point particle 
model (withp=l) which counts as a zero brane. Employment of brane models 
of lower dimension, p < n, for which the extrinsic confinement condition and 
the associated hyperbolicity requirement derived fom (7.2) are essential, is 
often useful for providing an approximate descriptions of higher dimensional 
case when the the fields characterising the latter are highly concentrated 
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in the neighbourhood of a lower dimensional world sheet within a distance 
that is small compared with the scales characteristic of dynamic variations 
in directions tangential to the world sheet. Thus for example a point particle 
model might be useful for describing the motion, with respect to a relatively 
slowly varying background, of a small loop in a string model that might itself 
be just a opproximation for decribing what at a more microscopically accu- 
rate level might need the use of a continuum model. The example that has 
been most important in motivating the development of the relativistic for- 
malism described here is that of the representation ( as originally suggested 
by Kibble'^^^1, Wittent^^^l and others) of vortex defects (due to spontaneous 
symmetry breaking) of the vaccuum by ( "cosmic" ) string models as a macro- 
scopic approximation for use in the (cosmologically important) cases in which 
the vortex thickness can be treated as negligible compared with other rel- 
evant length scales. This lead to the introduction of models of variational 
type in which the action was to be thought of as being derived from the 
microscopic action of the relevant underlying field theory by integral across 
the vortex in a local equilibrium state. 

Quite generally, in cases where a compound system has a variational 
formulation in terms of a total action of the form J2 ^, the action contribution 
of an individual p brane of the system will be given by a corresponding p 
surface integral 



where dS_ denotes the induced surface measure and £ is a Lagrangian scalar 
function of whatever internal fields on the world sheet are involved and also 
of any relevant externally induced fields such as those given by (7.2). In 
the simplest (non conducting) cosmic string models originally envisaged by 
Kibble'^^^l it was sufficient to use a Goto-Nambu'^^^^ action in which (as in 
the analogous Dirac membrane model^^^^') the scalar L is specified trivially as 
a constant, which is interpretable as the negative of the (spacially isotropic) 
tension T which in this case is not only uniform but (as an expression of 
the special property of intrinsic Lorentz invariance which distinguishes these 
particular models) is also equal in this case to the value of the energy density 
U say. The quantity U can be defined, for a generic brane model, as the 
eigenvectvalue corresonding to the timelike eigenvctor of its surface stress 
momentum energy tensor T^*^, while in a string model the tension T is un- 
ambiguously definable as the other eigenvalue, the case of a membrane being 
more complicated in far as it admits the possibility of two possibly distinct 




(7.3) 
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tension eigenvalues. 

It typically occurs that the approximate macroscopic treatment of a 
system that is conservative, with a variational formulation, at a microscopic 
level may require the use of a non conservative macroscopic model involving 
averaging over microscopic degrees of freedom that are taken into count as 
entropy. Although it may invalidate the conservative nature of the model as 
a whole, such an averaging process does not invalidate the local conservation 
laws obeyed by additive quantities such as energy momentum or electromag- 
netic charge : what happens is that instead of having the status of Noether 
identities expressing the invariance properties that hold for the underlying 
variational model, such conservation laws are to be interpreted in the macro- 
scopic model as constistency conditions for the existence of a corresponding 
microscopic variational model. The commonly but (not always) appropriate 
notion that a macroscopic model under consideration is obtainable by inte- 
grating out the fine details of a more complicated underlying model makes it 
seem physically natural to try to preserve some of the spirit of the original 
finer model by using a description in terms of Dirac distributions. However 
although very useful for some purposes when used with discretion, use of 
Dirac distributions can easily become addictive, and is often systematically 
abused in a manner that hinders clear analysis and provides an archetyp- 
ical example of the kind of excess mathematical baggage that the present 
approach is designed to avoid. 

The most important example of mathematical machinery that is very 
helpful for many specific purpose but whose use needs to be avoided (as 
excess baggage) when one wants to obtain a simple formulation of gen- 
eral principles such as that embodied in the "generalised sail equation" '^^1 
(1.1), is that of the introduction of a system of internal coordinates a* say, 
{i = 0, ...,p — l), on the p dimensional world sheet of the (p — l) brane under 
consideration, whose imbedding is thereby describable as a mapping a* i— > x'^ 
where the x^, (/x = 0, 1, n — 1) are local coordinates on the n dimensional 
background spacetime. A possibility that is of considerable practical util- 
ity in the intermediate stages of many calculation is the use of what I call 
"adapted coordinates" meaning a matched system of internal and external 
coordinates in terms of which the imbedding mapping is characterised by 
x" = a°, xP'^ = a^'^ , xp = 0, x"'^ = , but this obviously can 
not be done for the simultaneous treatment of intersecting branes (as at the 
junctions in a cluster of soap bubbles) and it is also obviously incompati- 
ble with the freedom to use an objective characterisation of the background 
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coordinates (e.g. in flat space applications by the requirement that they be 
Minkowskian) which may be important for the flnal presentation and utihs- 
abihty of the resuhs. 

One of the uses, as an intermediate step, of a coordinate mapping 
^ jg fQj. ^Yie exphcit construction of the corresponding intrinsic compo- 
nents of the images induced in the imbedding of covariant tensor flelds on the 
background space, such as the electromagnetic potential and most impor- 
tant of all the background space time metric g^i,, whose respective images 
are given by 

dx'^ , dx^ dx^ ,_ ^. 

A^^a, = A^— , g^,^h,,=g^,-^— . (7.4) 

In cases where a compound system has a variational formulation as a sum in 
which each distinct brane contributes a term of the form (7.3), the obvious 
analogue of the traditional variational speciflcation of the conserved current 
and stress energy momentum tensor (whose local conservation equations are 
the Noether identities expressing gauge invariance and general diffeomor- 
phism covariance) will take the form 

dC dC 

subject to the proviso (which is not necessary for the simple and conducting 
cosmic string models models that will be considered below)that Eulerian 
variational derivatives are to be used instead of simple partial derivatives 
if derivatives of the potential and metric are involved. The quantities h^^ 
appearing in (7.5) are of course the components of the contravariant inverse 
of the induced metric which is to be used for raising and lowering internal 
indices. 

Whether they are specifled variationally, as in (7.5), or whether they 
are specifled in some more empirical way, as would be necessary in a general, 
non conservative model, the internal current and stress energy momentum 
tensor will determine corresponding background tensor flelds by the natural 
pull back mapping that is determined directly by the imbedding for any 
contravariant vector flelds, the corresponding coordinate expressions being 
given by 
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The idea of the strategy developed here is that it is more efficient 
for general theoretical (as opposed to specific computational) purposes not 
to work with internal tensors such as and but rather to work with 
the corrsponding background spacetime tensors, which in this case are Jj^ 
(with the underlining as a reminder that it refersto a surface not volume 
current) and T^*^. When a variational specification is available it is preferable 
(particularly for dealing with compound systems involving several mutually 
interacting branes of diverse dimensions) to bypass the passage via (7.4) 
and (7.6) through the internal coordinate versions by replacing (7.5) by the 
equivalent but more direct background coordinate specifications 

in which the only formal difference from the usual expression for a continuous 
medium as opposed to lower dimensional brane model is the replacment in 
the last term of the contravariant version g^^ of the ordinary background 
metric by what I call the (first) fundamental tensor g^^'^ of the brane world 
sheet, which is obtained here as the pull back of the contravariant inverse of 
the induced metric, i.e. 

= , (7,8) 

where a double overline is introduced here to denote the surface tangential 
part of any tensor as defined with respect to the background metric, i.e. the 
result of contracting all its indices with respect to the mixed rankp projection 
tensor version ^^j, of the fundamental tensor itself. It is to be noted that 
any background tensor that is obtained as the puUback of an intrinsic tensor 
within the imbedded surface will automatically be equal to its own tangential 
part, so that in particular we shall have = J_^ and T_^'^ = T_^'^ . 

The "fundamental" tensor of the imbedding that is thus specified in 
accordance with (7.8) is of great (but still insufficiently widely recognised) im- 
portance as the starting point for the systematic tensorial analysis of imbed- 
ding curvature as described in the next section, whose results are applicable 
not just to a timelike brane world sheet but also to submanifolds that are 
spacelike (though not to those that are null, i.e. metrically degenerate). 

In the particular case of a Goto Nambu string model'^^^^ or a Dirac 
membrane model^^^^^, as characterised by an action of the form (7.3) with 
L = Lq for some fixed value Lq, which (see section 3) gives a uniform isotropic 
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tension T that is equal to the corresponding energy density U and opposite to 
the Lagrangian itself, i.e. U = T = —L^, the introduction of the fundamental 
tensor g^'^ makes it easy to check the well known property that the charac- 
teristic propagation speed of extrinsic perturbations is, in this case, that of 
light (c=l in the units used here) by substituting in (7.1) the simple formula 
whereby the stress momentum energy density for such a (Gotu Nambu or 
Dirac) model is expressible directly as Tj^^ = —Ug^^. 

As explained in the appendix, we shall adopt the systematic use of a 
convention using an overhead parallelism symbol, ^, to indicate the effect of 
projection into the surface, and an overhead perpendicularity symbol, — , to 
indicate the effect of the complementary orthogonal projection operation, so 
that the surface tangentiality conditions that the surface current and stress 
momentum energy tensors must satisfy by construction, will simply take the 
form 

= , = , (7.9) 

it can be seen from (7.7) that the variations in a brane Lagrangian L 
due to an infinitesimal electromagnetic gauge variation i— > + V^x ^md 
an infinitesimal diffeorphism variation g^^, g^y + V(^^,y) of the metric will 
be expressible respectively as 

J^V^X = ^,{xin - X^,I^ , (7.10) 

and 

r"V^6 = V^(6r") - e.V^r'^ . (7.11) 

It can be seen that the first term on the right of each of these equations 
has the form'^^1 that characterises a tangential current divergence within 
the p dimensional brane world sheet, and hence that by the appropriate p 
dimensional version of Green's theorem the corresponding surface integral 
will be expressible as the integral over the brane boundary (if any) of the 
contraction of the tangential current with the unit world sheet tangent vector 
normal to, and oriented towards, the boundary. 

Let us consider the very large class of situationst^^^ that can be rep- 
resented by a well behaved brane complex ( or "rigging system") in which 
direct action of a lower on a higher dimensional brane occurs only when the 
former forms a smooth boundary segment of the latter (as when a monopole, 
treated as a point particle, forms the termination of a string, or when a sail 
forms the boundary between two external wind volumes), subject to dynamic 
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equations to the effect that the infinitesimal variation of the relevant fields 
other than the externally determined background fields g^j^^, and A^, gives no 
contribution to the variation of the combined action J2 ^ taken over the var- 
ious brane constituents of the system, restricting ourselves to cases in which 
derivatives of the external fields gfj_i, and are not involved in the action. 
(The exclusion of more general derivative couplings merely avoids the extra 
technical complications that are present in more elaborate, e.g. polarised sys- 
tems, but the exclusion of direct action except on a smooth boundary is more 
essential, being needed to avoid the serious divergence difficulties, exempli- 
fied by that of the radiation back reaction on a point particle, which would 
otherwise be involved.) Then it can be seen (by systematically using (7.10) 
to convert divergences to boundary contributions) that the requirement that 
this combined action should also be identically invariant under gauge 
transformations generated by an arbitrary field x is equivalent to the condi- 
tion that there should be a total current conservation law expressed by the 
condition'^^^l that for each p brane of the system we should have 

where the summation is taken over the separate (p+ 1) branes of which the p 
brane under consideration forms part of the bondary, and where J'^ without 
underline denotes the value on the boundary segment of the current vector 
in the higher dimensional sheet while denotes the unit normal from the 
p dimensional boundary into the relevant externally attatched brane world 
sheet. Similarly (by analogous systematic use of (7.11) to convert divergences 
to boundary contributions) it can be seen under the same conditions that the 
general covariance requirement that the combined action be invariant under 
diffeomorphisms generated by an arbitrary vector field is equivalent to a 
local energy momentum conservtion law to the effect that for each brane 
of the system we should have 

V.T'^^ = r, f, = Y.^'^T\ + F,,r (7.13) 

in which the force density is obtained as the sum of contact contributions 
from the (non underlined) stress momentum energy density tensor T^^'^ of 
each of the attached p branes (at most two if p = n, but arbitrarily many for 
p < n) oi which the (p — 1) brane under consideration is a boundary segment, 
together with an external electromagnetic force contribution determined by 
the Maxwelhan field F^j/ = 2V[^A^]. 
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Although the foregoing direct derivation starts from a variational pos- 
tulate, charge and energy momentum conservation laws of the form (7.12) 
and (7.13) can still be expected to hold for more general dissipative models 
such as would be obtained by macroscopic averaging over internal degrees of 
freedom whose net effect would be taken into account in terms of entropy 
currents. An alternative (for some tastes more intuitive, but mathematically 
much more awkward) way of deriving (3.4) and (3.5) in such cases would 
be to consider the brane system as the infinitely thin limit of a continuous 
medium model where the current and stress energy momentum density 
T^^'^ are no longer continuous fields but have become Dirac distributions, 
whose coefficients are interpretable as the corresponding smooth world sheet 
supported fields J_^ and Tj^^ . By whatever route they may have been ob- 
tained, the ubiquitous generality of (7.12) and (7.13) - and of the extrinsic 
equation of motion (7.1) that is obtainable via (7.9) as a direct consequence 
- cannot be overemphasised. In the particular case of a free motion for which 
external electromagnetic and contact effects are absent we evidently get 

^ = ^ T^''K^/ = {) (7.14) 

In the case of a variational model with action simply proportional to the 
world sheet measure, as in the case of a Dirac membrane, a Goto-Nambu 
string, or an ordinary free point particle, the force free equation of extrinsic 
("brane wave") motion (7.14) obviously reduces to the even simpler ("har- 
monic") form K^^ = which includes the equation for a geodesic in the one 
dimensional case. 

8. Perfect Brane Models. 

For a general brane model, we can always define an energy density 
scalar, U say, as the negative of the eigenvalue specified by 

T'^^u'' = -Uu" (8.1) 

where the corresponding eigenvector u'^ is distinguished by the requirement 
that it be timelike or null. As a widely applicable special case (including the 
Dirac membrane mentioned above, as well as all point particle and string 
models) a {p — 1) brane may be described as "perfect" if its surface stress 
momentum energy tensor is isotropic with respect to the other orthogonal 
directions, which in the generic case for which the eigenvector u'^ is strictly 
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timelike (not null) and hence normalisable to unity, one gets'^' the explicit 
form 

T^, = {U -T)u^u,-TT, , n% = -l, (8.2) 

where T (the negative of the other (p—l) degenerate eigenvalues) is what is 
interpretable as the tension of the {p — 1) brane. 

The category of perfect branes includes, as the extreme case p = n, 
the example of an ordinary "perfect fluid" (with U = p, where p is the ordi- 
nary volume density of mass-energy, while T = —P where P is the ordinary, 
positive, pressure). In the other cases, i.e. for a {p — 1) brane of lower dimen- 
sion than the background, i.e. p < n, for which extrinsic displacements are 
possible (so that the tension must be non negative in order to avoid local 
instability'^^1'^^^1) the extrinsic motion will be governed by (7.1) or in the 
force free case by (7.14) which, on substitution of (8.2) gives the dynamic 
equations for a free perfect brane world sheet in the form 

cjK'^ = (1 - c/)i^n'^ , = u'^V.u^ , = (8.3) 

where ii^ is the acceleration vector of the unit eigenvector and is inter- 
pretable as the speed of propagation - relative to the preferred frame specified 
by - of extrinsic perturbations, as derived from the general characteristic 
equation (7.2). It is to be noted that in the ultra relativistic case of a Dirac 
membrane or Goto Nambu string one has = 1 which means that the 
right hand side of (8.3) will vanish. On the other hand the strings and mem- 
branes that are commonly used (in violins, drums, etc.) by old fashionned 
non relativistic (i.e. non electronic) orchestras for music generation, will also 
be describable to a very good approximation by this same equation but with 
<< 1, which means that the coefficient c^^ will be able to be neglected 
on the right though not of course on the left. 

The extreme case of a "zero brane" with p = 1, i.e. that of an 
ordinary (massive) point particle, can be considered as being automatically 
of the perfect type characterised by (8.2) with U = m where m is its mass, and 
with identically vanishing tension T = which is consistent with the obvious 
necessity of having zero relative speed of propagation of any perturbation in 
this one dimesional case. For a point particle trajectory the first and second 
fundamental tensors will be given simply by 

^/^^ = -u>'u, , K^/ = ^^,KP , = -M^ , = u^y.u^ (8.4) 
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while in terms of the particle mass m and charge e say substitution of the 
appropriate expressions 

T^^ = -m , = eu^ , (8.5) 

into the general expressions (7.12) and (7.13) gives the dynamical equations 
in the familiar form 

M^V^e = , u^V^m = , -mK^ = eF^.u^' , (8.6) 

subject of course to the usual proviso (which in this context is to be taken 
quite litterally!) that there are no strings attached, since otherwise corre- 
sponding contact contributions on the right of (7.12) and (7.13) could cause 
variations of the values of the charge and mass scalars, e and m as well as 
modifying the acceleration equation in (7.14). 

The case of a membrane in 4-dimensions (or more generally of 
an (n — 2) brane in n dimensions) shares with the opposite extreme case 
of a point particle the property of having comparatively simple kinematic 
properties, since any timelike hypersurface has first and second fundamental 
tensors that are expressible in terms of its unit normal (as specified by an 
arbitrary choice of orientation) in the form 

Tu = g\ - , K^/ = K^^y , A^A^ = 1 . (8.7) 

Analogously to the way the first fundamental tensor 'g^'^ is specifiable (by 
(7.8)) as the pull back of the contravariant version of the induced metric, i.e. 
of what is commonly known as the first fundamental form of the imbedding, 
so analogously the symmetric tensor K^^ is the pull back of the contravariant 
version of what is commonly known as the second fundamental form on the 
hypersurface, a quantity whose specification, like that of the unit normal A^ 
involves an arbitrary choice of sign. (In addition to its principle advantage 
of being applicable to imbeddings of arbitrary dimension, not just hypersur- 
faces, an advantage of our present strategy of working with the three index 
second fundamental tensor rather that the two index second fundamental 
form even in the hypersurface case where the latter is available is that unlike 
that of K^u the specification of K^^p is quite unambiguous.) Whereas the 
kinematic specifications (8.7) are simpler than their analogues for the lower 
dimensional case of a string, on the other hand the dynamics of a membrane 
are generally more complicated. Unlike the case of a string model which must 
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always, trivially, be perfect in the sense of (8.2) (or of its null limit t^^l) the 
postulate of "perfection" in this sense is a serious restriction in the case of a 
membrane, being satisfied for a Dirac membrane or an ordinary soap bubble 
type membrane, (and even as a reasonable approximation to the way musical 
drum membranes are most commonly tuned), but it will not be at all valid 
for such applications as to a typical ship's sail. 

Between the highpersurface supported case of a membrane and 
the curve supported case of a point particle the only intermediate kind of 
brane that can exist in 4-dimensions is that of 1-brane, i.e. a string model, 
which (for any background dimension n) will have a first fundamental tensor 
that is expressible as the square of the antisymmetric tangential tensor S'^'^ 
that is defineable'^^^1 as the puUback of the contravariant version of the in- 
duced measure tensor that is specified modulo a choice of orientation by the 
imbedding, i.e. we shall have 

= S^S'u , = ^f^"' . (8.8) 

A special feature distinguishing string models from point particle 
models on one hand and from higher dimensional brane models on the other 
is the dual symmetry t^^^"^^' that exists at a formal level between the spacelike 
and timelike eigenvectors u'^ (as already introduced) and v'^ that for a generic 
case (excluding the null state limitt^^^) are characterised modulo a choice of 
orientation by the expression 

X!"" = Uu^u" - Tvf'v" , yf'Vf, = 1 = -uf'u^ (8.9) 

in which the tension T appears as the dual analogue of the "rest frame" 
energy per unit length U. This formal duality can also be made apparent in 
the expression for the extrinsic curvature vector of the string, which can be 
expressed as 

K'' = -^"uiv" - ii") , v'^" = v'Vuv" , ii'' = v"^,u^ (8.10) 

whose substitution in (8.3) enables the equation of extrinsic motion of a free 
string to be expressed in the manifestly self dual form 

U^^ii" = t\\v'' . (8.11) 

Of course the extrinsic equation of motion, whether of the general 
form (7.1) or the free string specialised form (8.11), cannot actually be used to 
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determine the evolution of the world sheet until the appropriate prescription 
has been given for evaluating the necessary stress momentum energy tensor 
components, which in the string case (8.11) can be taken to be just T and 
U. In the simple Goto-Nambu case, for which these eigenvalues are specified 
in advance to have constant values, U = T = —L^, no further preparation 
is needed for the integration of (8.12) but in general, for a string model 
with non trivial intrinsic structure the completion of the system of equations 
of motion will involve the specification of other differential equations. The 
simplest non trivial possibility, which is applicable to higher dimensionsional 
perfect brane models as well as to strings, is what is known in the specific 
context of perfect fluid theory as the "barotropic" case, meaning the case in 
which T is specified (directly or parametrically) as a function only of [/ by a 
single equation of state. In this barytropic case (which includes the Witten 
type conducting cosmic string modelst'^^^' whose investigation provided the 
original motivation for this work) the only differential equations that are 
needed to supplement the extrinsic equation of motion (8.3) or (8.11) are 
those that are obtained from the projection into the world sheet of the full 
local momentum energy conservation equation (7.13), which in the force free 
case simply gives 

W^^ = (8.12) 

whose two independent components can be conveniently expressed as a pair 
of mutually dual surface current conservation laws given by 

V^(i^n^) = , V^ifLv'') = , (8.13) 

in terms of an effective number density u and an associated effective mass 
density fi that obtained from the equation of state as functions of U or 
equivalently of T by a pair of (mutually dual) integral relations of the form 

In = / , In u = / , (8.14) 

J U-T ' ^ J T-U' ^ ' 

which fix them modulo a pair of constants of integration of which one is 
conventionally fixed by imposing the (self dual) restraint condition 

^u = U-T. (8.15) 

Appart from the extrinsic perturbations of the world sheet loca- 
tion itself, which propagate with the "brane wave" speed ce (relative to the 
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frame deterined by m'^) as afready discused, the only other kind of pertur- 
bation mode that can occur in a barytropic string are longitudinal modes 
specified by the variation of U or equivalently of T within the world sheet. 
Such longitudinal perturbations (the analogue of ordinary sound waves in a 
perfect fiuid) can easly be seent^^lt^^^l to have a relative propagation velocity 
given by 

which must be real in order for the string to be locally stable. Knowledge 
of whether the longitudinal perturbation speed is greater or less than the 
extrinsic speed may be critically significant for questions such as the sta- 
bility of stationary rotating ring equilibrium states'^^lt^^^ll^^^^l^^^' and their 
deformed generalisations'^^^' which for Witten type cosmic strings (as op- 
posed to the ordinary Goto Nambu type for which no such states exist) 
may be cosmologically important '^^^"^^^1. Most early, and many more recent 
discussions[^^^l[i^5][^^6l[i^^] of Witten type strings were implicitly based on the 
use of an equation of state for which the sum U + T remains constant, which 
implies longitudinal propagation at a speed equal to that of light, = 1 
which thus necessarily exceeds but more accurate investigations'^^^^t^^^' 
have recently been developed'^^^^ to a stage at which it is becoming increas- 
ingly clear that the opposite is usually the case, i.e. Witten type models 
would seem to be typified by < c^. 

A very special interest attaches to the intermediate non-dispersive, 
case characterised by = , which corresponds to an equation of state for 
which the eingenvalue product TU is constant, leading to dynamic equations'^^^' 
that I have shown to be explicitly integrable (like those of the degenerate Goto 
Nambu case) in a flat spacetime background, the general form in an arbitrary 
curved background being expressible as 

where is a constant and the (timelike) unit vectors L±'^ are directed along 
the "left" and "right" moving unit characteristic directions, the former be- 
ing parallel propagated by the latter and vice versa. Another, more recently 
established special property of the non-dispersive equation of state UT = 
is the remarkable solubility (by separation of the relevant Hamilton Jacobi 
equation) of the corresponding equations of stationary equilibrium not only 



82 



in flat space but in a generalised Kerr black hole background'^^"' . This spe- 
cial "constant product" string model (which can be recognised '^^^ as turning 
up spontaneously in Kaluza Klein theoryt^^^lt^^^l'^^^^t^^^l is not just of purely 
mathematical interest : my predctiont^^^l that it should provide a good de- 
scription of the averaged effect of random noise perturbations on an "ordi- 
nary" Goto-Nambu type cosmic string (on the grounds that their presence 
should not introduce dispersion) has been confirmed by Vilenkin's more de- 
tailed "wiggly string" calculations^^^^^ . 

Appendix : Background tensor analysis of the curvature of an 
imbedding. 

In any n dimensional manifold with a non degenerate (Rieman- 
nian or pseudo Riemannian) metric tensor with components g^j^y (with re- 
spect to some local coordinate patch) that is to be used for index lowering 
and raising, any non-null (strictly spacelike or timelike) p-dimentional surface 
element at a point determines a corresponding decomposition 

where 'g^y is the fundamental (rank p) projection tensor of the surface ele- 
ment, and ~g^y is the complementary (rank n — p) tensor of projection or- 
thogonal to the surface. Consistently with (Al) we shall adopt the systematic 
use of a convention using an overhead parallelism symbol, ^, to indicate the 
effect of projection into the surface, and an overhead perpendicularity sym- 
bol, — , to indicate the effect of the complementary orthogonal projection 
operation, so that in particular for an arbitrary vector with components 
and for the standard operator of Riemannian covariant differentiation (as 
defined with respect to a symmetric connection such that V^g^p vanishes), 
we have 

def —— def I def i I def i 

e = TuC , = r^v. , 1^ = , = t%v, . {A2) 

In terms of this convention, the fundamental tangential and orthogonal pro- 
jection operators are thus characterised by the conditions that for any vector 
that is tangent to the p-surface element, and any vector A'* that is orthog- 
onal to the p-surface element we must have 

W = u^ , -b^ = , = , \^ = V . (A3) 
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Unlike the full covariant differentiation operator and its or- 
_i_ 

thogonally projected part V^, the tangential covariant differentiation oper- 
ator, V^, has the property of being well defined not only for (sufficiently 
smooth) fields defined on an open neighbourhood of the background space 
but even for fields with support is confined to a (sufficiently smooth) p-surface 
whose tangent surface element specifies the projection. In particular, for any 
such p-surface there will be a well defined second fundamental tensor, K^^'^ 
definedl^^^ll^^l'^^l in terms of its first fundamental tensor ^^j, by 

which as a trivial algebraic identity is obviously tangential on the first two 
indices and almost as obviously orthogonal on the last, i.e. for an arbitrary 
vector C,^ it satisfies 

K^/^, = = K^Ap ■ (Ab) 

Such a tensor K^iy'^ is of course definable not only for the funda- 
mental projection tensor of a p-surface, but also for any (smooth) field of 
rank p projection operators g^^i, as specified by a field of arbitrarily orien- 
tated p-surface elements. What distinguishes the integrable case, i.e. that in 
which the elements mesh together to form a well defined p-surface through 
the point under consideration, is the condition that the tensor defined by 
(A5) should also satisfy the Weingarten identity 

(where the square brackets denote antisymmetrisation), this non trivial sym- 
metry property of the second fundamental tensor being derivable '^^^1 as a 
version of the well known Frobenius theorem. 

The second fundamental tensor K^^^ has the property of fully 
determining the tangential derivatives of the first fundamental tensor ^^j, by 

the formula 

V,J.p = 2K,,(,p) (AT) 

(using round brackets to denote symmetrisation) and it can be seen to be 
characterisable by the condition that the orthogonal projection of the accel- 
eration of any tangential vector field will be given by 

-^^u'V^u^ = u^u-'K^/ , (AS) 
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as well as by the condition (in which the non-trivial role of the symmetry 
property (A6) is more apparent) that the tangential projection of the deriva- 
tive of any field of surface normal vectors should be given by 



V^A, = -i^^/A, . {A9) 

Going on to higher order we can introduce the third fundamental 
tensor'^^] in an analagous manner as 

Sam/ = T^.TuTu-l'aVxK^r" , (AlO) 

which by construction is obviously symmetric between the second and third 
indices and tangential on all the first three indices while being, i.e. (for an 
arbitrary vector C,'^) it satisfies the trivial identities 



In a spacetime background that is flat (or of constant curvature as is the 
case for the DeSitter universe model) this third fundamental tensor is fully 
symmetric over all the first three indices by what is interpretable as the 
generalised Codazzi identity which is expressible t^^^^ in a background with 
arbitrary Riemann curvature R\p^a as 

Sam/ = ^{Xpuf + lWxW{pTu)R.r^erg'(, • (^12) 

It is very useful for a great many purposes to introduce the ex- 
trinsic curvature vector K^, defined as the trace of the second fundamental 
tensor, i.e. 

= KY . T^' = {A13) 

The specification of this extrinsic curvature vector for a timelike p-surface in 
a dynamic theory provides what can be taken as the equations of extrinsic 
motion of the p-surface'^^^ (the simplest case being the "harmonic" condition 
K'^ = obtained from a simple surface measure variational principle such 
as that of the Goto-Nambu string model or the Dirac membrane model). It 
is also useful for many purposes^^^' to introduce the extrinsic conformation 
tensor C^^'^ defined as the trace free part of the second fundamental tensor 
by 

def 
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which (hke the Wey tensor of the background metric) has the noteworthy 
property of being conformally invariant with respect to conformal modifica- 
tions of Qf^u ^ &'^"9iJiv of the background metric. 

The condition of preserving the tangent element to an imbed- 
ded p-surface at a point breaks down the full n dimensional rotation group 
preserving the background metric into the product of the restricted p dimen- 
sional rotation group preserving the induced metric in the imbedding with 
the restricted (n—p) dimensional rotation group preserving the induced met- 
ric in the orthogonal element. Associated with each of these subgroups there 
is a corresponding naturally induced connection and covariant differentiation 
operator acting on the corresponding bundles of tangent vectors and or- 
thogonal vectors A'^ respectively, and for each there will be a corresponding, 
respectively "inner" and "outer" bundle curvature, which will be represented 
by a corresponding background tensor, the former "inner" curvature tensor 
being just the pull-back onto the background by the imbedding mapping of 
the ordinary Riemann curvature of the intrinsic geometry induced by the 
imbedding. Explicitly for any vector fields satisfying the appropriate tangen- 
tiallity and orthogonality conditions (A2), the effects of the corresponding 
restricted "inner" (tangentially projected) and "outer" (orthogonally pro- 
jected) differentiation operations will be given respectively by 



V^f^ = Tp^f.u' = - K^p^u" , TpV^A" = V^A'^ + K^pY . {Alb) 

Using the convention that an underline is inserted whenever necessary to 
distinguish a quantity defined with respect to an imbedding from its higher 
dimensional background analogue, the corresponding inner curvature tensor 
RkX^v of the p-surface (as distinct from the ordinary background Riemann 
tensor R^x^u) and the corresponding outer curvature tensor Qkx'^u (for which 
no background analogue exists, so that there is no need to underline it) are 
specifiable by the respective conditions that for any tangential vector and 
any orthogonal vector A^' to the surface, i.e. for any vectors satisfying (A3) 
we should respectively have 



2V[^V,]nP = 2^^r[.V^](rpV,n,) = Ep/<x«" , (^16) 
and 

2VAr[.VH(i%v,A,) = ^]^/.A^. (Air) 



86 



Then it can be verified that the inner curvature tensor is given in terms of 
the tangential projection of its background analogue by the relation 

R^,P„ = 2A'^[/A',]., + 1^/^ , (A18) 

which is the translation into the present scheme of what is well known in other 
schemes as the generalised Gauss identity. The much less well known analogue 
for the (identically trace free and conformally invariant) outer curvature, for 
which the most historically appropriate name is arguably that of Schouten, 
is given^^^l by the expression 

= 2Cy^^'C,^r. + r#.aA%i''«V<x , (^19) 

where 0^^^^ is (trace free conformally invariant) background Wey tensor, 
which is definable implicitly for a background of dimension n > 2 by the 
decomposition of the Riemann tensor into trace and trace free parts as 

where, as usual the background Ricci tensor and Ricci scalar are given by 

Rtiu = Rp/y , R = R^u ■ (^21) 

It can be seen from the form of the identity (A19) that in a flat or confor- 
mally fiat background (for which it is necessary, and for n > 4 sufficient, that 
the Weyl tensor should vanish) the vanishing of the extrinsic conformation 
tensor 0^,^^ will be sufficient (independently of the behaviour of the extrin- 
sic curvature vector K^) for vanishing of the outer curvature tensor VL^y^^, 
and hence (by (A17)) for the possibility of constructing fields of orthogonal 
vectors that satisfy the generalised Fermi- Walker propagation condition 

to the effect that ~g'' p^^Xp should vanish. It can also be shown (taking spe- 
cial trouble for the case p = 3 ) that in a conformally flat background (of 
arbitrary dimension n) the vanishing of the conformation tensor C^j^yP is al- 
ways sufficient (though by no means necessary) for conformal flatness of the 
induced geometry in the imbedding. 
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